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•The  study  is  part  of  the  DNA  C-6  Program  to  investigate  the  feasibility  of 
one  of  the  multiple  aim  point  (MAP)  basing  concepts  for  the  advanced  M-X 
missile  systems;  namely,  the  shelter  concept.  Its  objective  is  to  Identify 
requirements  in  the  finite  element  analysis  of  the  shelter  response  when 
the  structure  is  subjected  to  nuclear  environments,  and  to  define  an 
analysis  methodology  which  can  be  followed  in  v'erforming  dynamic  analyses 
of  a shelter-like  structure.  ^ 
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SUMMARY  AND  RECOMMENDATIONS 


This  study  examines  the  application  of  dynamic  finite  element  methods  to  shelter- 
like structures  In  a nuclear  environment.  In  particular,  the  study  investigates  the  deforma- 
tion modes  of  the  shelter  under  different  loading  conditions,  their  sensitivity  to  varia- 
tions in  the  overpressure  wave  form,  important  features  of  the  shelter  which  affect  the 
response,  the  effect  of  inelastic  material  properties  and  the  role  of  elastic  analyses  and 
ground  motions  in  a shelter  analysis  of  this  kind. 

Based  on  results  of  case  studies  obtained  using  implicit  elastic  (GENSAP)  and 
explicit  inelastic  (TRANAL)  finite  element  codes  an  analysis  methodology  is  refined  and  a 
procedure  for  performing  dynamic  analysis  of  the  shelter  is  developed.  The  procedure  in- 
volves modeling  the  shelter-like  structure  and  its  medium  according  to  the  guidelines  es- 
tablished in  the  study,  factoring  in  those  characteristics  of  the  air  overpressure  which 
have  been  found  significant  and  separating  the  analysis  into  two  distinct  but  complementary 
components;  short-term  analysis  and  long-term  analysis. 

In  the  short-term  analysis  the  response  of  the  front-end  of  the  shelter  which  is 
exposed  to  direct  airblast  load  is  examined  using  a refined  model  which  includes  detailed 
subcomponents  of  the  front  end  such  as  the  closure,  gap,  bearing  ring,  etc.  The  portion  of 
the  shelter  aft  of  the  transition  section  may  be  replaced  by  a simpler  model  which  approxi- 
mates its  impedance  and  energy  radiation  characteristics.  This  analysis  yields  peak 
stresses  in  the  headworks,  frame  and  closure  and  peak  deformation  of  the  closure. 

In  the  long-term  analysis,  a coarse  model  of  the  front  end  is  used  to  approxi- 
mate the  front  load  distribution  and  its  transfer  to  the  tubular  section.  The  long-term 
analysis  which  yields  the  gross  motion  of  the  headworks  and  the  tubular  section  is  best 
carried  out  using  the  soil  island  approach  to  minimize  boundary  effects. 

This  methodology  is  a natural  consequence  of  the  shelter  concept  wherein  the 
structure  is  subjected  both  to  direct  airblast  which  dominates  short-term  response  of  the 
front  end  and  indirect  ground  motion  which  is  responsible  for  the  long-term  gross  motion 
of  the  shelter.  The  distinction  between  short  and  long-term  response  is  also  dictated 
by  computational  expediancy;  this  combination  makes  comprehensive  three-dimensional  analyses 
of  the  shelter  structure  possible  and  economically  acceptable. 

The  present  study  contributes  to  the  definition  of  the  methodology  and  substan- 
tiates the  recommended  approach  with  results  from  case  studies.  Its  emphasis,  however, 
has  been  on  the  short-term  response  of  the  structure.  Even  there  the  scope  of  the  work  is 
limited  by  the  airblast  data  available.  Limited  attention  has  been  devoted  to  the  other 
facet  of  the  methodology,  l.e.,  the  long-term  response  and  especially  the  role  of  the  soil 
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Island  approach.  Observations  made  in  this  study  regarding  the  long-term  response  should  | 

be  confirmed  by  additional  studies  at  least  as  extensively  as  the  conclusions  on  short-  - 

term  response  have  been  substantiated.  1 
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element  analysis  of  the  shelter  response  when  the  structure  Is  subjected  to  nuclear 
environments  and  to  define  an  analysis  methodology  which  can  be  followed  In  performing 
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CONVERSION  FACTORS 
FOR  U.S.  CUSTOMARY  TO  METRIC  (SI) 
UNITS  OF  MEASUREMENT 


To  Convert  From  To  Multiply  By 


degree  (angle) 

radian  (rad) 

1.745 

329 

X E -2 

foot 

meter  (m) 

3.048 

000 

X E -1 

foot-pound-force 

joule  (J) 

1.355 

818 

Inch 

meter  (m) 

2.540 

000 

X E -2 

pound-force  (Ibf  avoirdupois) 

newton  (N) 

4.448 

222 

pound-force  inch 

newton-meter  (N*m) 

1.129 

848 

X E -1 

pound-force/inch 

newton/meter  (N/m) 

1.751 

268 

X E +2 

2 

pound-force/ foot 

kilo  pascal  (kPa) 

4.788 

026 

X E -2 

2 

pound-force/inch  (psi) 

kilo  pascal  (kPa) 

6.894 

757 

pound-mass  (Ibm  avoirdupois) 

kilogram  (kg) 

4.535 

924 

X E -1 

pound -mass-foot  (moment 

2 2 

of  inertia) 

kilogram-meter  (kg*m  ) 

4.214 

Oil 

X E -2 
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I 

II 


III 


TABLE  OF  CONTENTS 


INTRODUCTION  

SCOPE  

2.1.  PHASE  ONE  

2.2.  PHASE  TVJO 

PHASE  1 RESULTS  

3.1.  CALCULATION  lA— TWO-DIMENSIONAL  ELASTIC  CALCULATION  WITH 

SIDE-ON  INCIDENCE  

3.2.  CALCULATION  IB— TWO-DIMENSIONAL  ELASTIC  CALCULATION  WITH 

HEAD-ON  INCIDENCE  

3.3.  CALCULATION  IC— THREE-DIMENSIONAL  ELASTIC  CALCULATION  WITH 

HEAD-ON  INCIDENCE  AND  FRONT-FACE  GAP  

3.4.  CALCULATION  ID— THREE-DIMENSIONAL  ELASTIC  CALCULATION  WITH 

30°  OBLIQUE  FRONT-ON  INCIDENCE  


Page 

1-1 

2-1 

2-1 

2- 3 

3- 1 

3-1 

3-3 

3-6 

3-8 


IV  INTERPRETATION  OF  PHASE  1 RESULTS  4-1 

4.1.  DEFORMATION  MODE 4-1 

4.2.  COMPUTER  MODELING  CONSIDERATIONS  4-2 

4.3.  AIRBLAST  CHARACTERISTICS  AND  SHELTER  RESPONSE  4-3 

i 

V PHASE  2 RESULTS 5-1 

5.1.  CALCULATION  2A— THREE-DIMENSIONAL  INELASTIC  CALCULATION 

WITH  HEAD-ON  INCIDENCE  5-2 

5.2.  CALCULATION  2 B— THREE-DIMENSIONAL  CALCULATION  WITH  HEAD-ON 

INCIDENCE;  INELASTIC  MODEL  FOR  CONCRETE  ONLY  5-6 

5.3.  CALCULATION  2C— THREE-DIMENS lONAL  ELASTIC  CALCUIJVTION  WITH 

HEAD-ON  INCIDENCE  5-8 

VI  INTERPRETATION  OF  PHASE  2 RESULTS  6-1 

6.1.  DEFORMATION  MODE 6-1 

6.2.  STRESS  DISTRIBUTION  IN  SHELTER  6-3 

VII  CONCLUSIONS 7-1 

7.1.  STRUCTURAL  MODELING  CONSIDERATIONS  7_1 

7.1.1.  Closure  7_1 

7.1.2.  Frame  and  Gap 7-2 

7.1.3.  Headworks  7-2 

7.1.4.  Transition  and  Tube 7_2 

7.2.  DEFORMATION  MODE 7_2 


5 


r 

* 

I 


4 

I TABLE  OF  CONTENTS  (CONTINUED) 

i Section  Page 

! 

' 7.2.1.  Head-On  Loading  7-3 

7.2.2.  Oblique  Front-On  Loading  7-4 

7.2.3.  Side-On  Loading  7-4 

7.3.  AIRBLAST  CHARACTERISTICS  AND  DEFORMATION  MODES  7-4 

7.3.1.  Spatial  Characteristics  7-4 

7.3.2.  Temporal  Characteristics  7-5 

7.4.  EFFECT  OF  INELASTICITY 7-6 

7.4.1.  Elastic  Versus  Inelastic  Analysis  7-6 

7.5.  GROUND  MOTION  EFFECT  7-7 

7.6.  SHELTER  ANALYSIS  METHODOLOGY  7-8 

VIII  REFERENCES 8-1 

APPENDIX  A— MODIFICATIONS  TO  TRANAL  A-1 

A.l.  APPLIED  PRESSURE  LOAD  MODIFICATIONS  A-1 

A. 1.1.  Overall  Loading  A-1 

A. 1.2.  Detail  Loading  A-1 

A. 1.3.  Code  Modifications  A-2 

A. 1.4.  Sample  Check  Problem  A-4 

A. 2.  OTHER  MODIFICATIONS  A-4 

A.  2.1.  Region  Coordinate  Edge  Processing A-4 

A. 2. 2.  Stability  and  Volume  Check  A-5 

A. 2. 3.  Region  Material  Description  A-5 


LIST  OF  FIGURES 


3-12b 

3-13 

3-14 

3-15a 

3- 15b 

3-16 

3-17 

3-18 


SI  shelter  configuration  

Two-dimensional  (plane  strain)  finite  element  model  of  shelter 
and  media,  slde-on  incidence.  Calculation  lA  

Comparison  of  peak  overpressure  profiles,  two-dimensional  slde-on 
incidence.  Calculation  lA  

Comparison  of  pressure/time  histories,  free-fleld,  slde-on  incidence 
Calculation  lA  

Comparison  of  pressure/time  histories,  berm,  side-on  Incidence, 
Calculation  lA  

Deformation  mode  T = 40  msec.,  side-on  incidence.  Calculation  lA  . . 

Deformation  mode  at  T = 60  msec. , slde-on  Incidence,  Calculation  lA 

Fourier  spectrum  of  deformational  mode  at  T = 40  msec. , side-on 
incidence.  Calculation  lA  

Fourier  spectrum  of  deformational  mode  at  T = 60  msec.,  side-on 
incidence.  Calculation  lA  

Hoop  compression  and  bending  stress  time  histories 

crown,  slde-on  Incidence,  Calculation  lA  

Soil  response,  top  of  berm,  side-on  Incidence,  Calculation  lA  . . . 

Horizontal  and  vertical  stress/time  histories  in  soil,  GZ  spring- 
line, side-on  incidence.  Calculation  lA  

Horizontal  and  vertical  stress/time  histories  in  soil,  rear 
springline,  slde-on  Incidence,  Calculation  lA  

Three-dimensional  finite  element  model  of  shelter  and  media, 
front-on  Incidence,  Calculations  IB  and  1C  

Cross-sectional  view  of  three-dimensional  model,  front  face 

and  headworks  

Cross-sectional  view  of  three-dimensional  model,  transition 
section  and  tube  

View  of  the  three-dimensional  model  of  the  bare  shelter  structure, 
soil  elements  not  shown  

Comparison  of  peak  overpressure  profiles,  front-on  incidence. 
Calculations  IB  and  1C  

Comparison  of  overpressure  waveforms,  near  front  face,  at  front 
face,  front-on  incidence  

Comparison  of  overpressure  waveforms,  top,  front-on  incidence  . . . 

Comparison  of  peak  overpressure  distribution  on  front  face  

Deformation  modes,  exaggerated,  AFWL  8041  load,  front-on 
incidence.  Calculation  IB  

Axial  stress-time  histories  along  shelter  length,  upper  head- 
works,  front-on  incidence.  Calculation  IB  


LIST  OF  FIGURES  (CONTINUED) 


Figure  No.  Page 

3-19  Axial  stress/time  histories  along  shelter  length,  lower  headworks, 

front-on  incidence.  Calculation  IB  3-33 

3-20  Axial  stress  time  histories,  tubular  portion  behind  headworks, 

front-on  Incidence,  Calculation  IB 3-34 

3-21  Axial  stress  time  histories,  tubular  portion  half  tube  diameter 

behind  headworks,  front-on  incidence.  Calculation  IB  3-35 

3-22  Distribution  of  peak  bearing  load  exerted  by  closure  on  headworks, 

front-on  incidence.  Calculation  IB  3-36 

3-23  Deformation  mode,  exaggerated,  gap,  front-on  incidence.  Calculation  1C  ...  . 3-37 

3-24  Deformation  pattern  of  closure  and  front  face  at  t = 10  msec.  , 

exaggerated,  front-on  incidence.  Calculation  1C  3-38 

3-25a  Plot  of  first  and  second  stress  invariants  at  front  face,  gap, 

front-on  incidence.  Calculation  1C  3-39 

3-25b  Plot  of  first  and  second  stress  invariants  at  front  face,  gap, 

front-on  incidence.  Calculation  1C  3-40 

3-26  Effect  of  gap  on  peak  bearing  load  exerted  by  closure  on  headworks 3-41 

3-27  A view  of  finite  element  model  for  quartering  calculation  (30°) , 

Calculation  ID 3-42 

3-28  Deformation  mode,  vertical  plane  of  symmetry,  30°  oblique 

front-on  incidence.  Calculation  ID  3-43 

3-29  Deformation  mode,  horizontal  plane  of  symmetry,  30°  oblique 

front-on  incidence.  Calculation  ID  3-44 

3-30  Variation  in  axial  stress  along  top  of  structure,  30° 

oblique  front-on  Incidence,  Calculation  ID  3-45 

3-31  Variation  in  axial  stress  along  bottom  of  structure,  30° 

oblique  front-on  incidence.  Calculation  ID  3-46 

3-32  Variation  in  axial  stress  along  springline,  blast  side 

of  structure,  30°  oblique  front-on  incidence.  Calculation  ID  3-47 

3-33  Variation  in  axial  stress  along  springline,  back  side  of 

structure,  30°  oblique  front-on  incidence,  Calculation  lA  3-48 

3-34  Axial  stresses  at  transition,  difference  in  axial  forces 

reflects  bending,  30°  oblique  front-on  incidence.  Calculation  ID  3-49 

3- 35  Comparison  of  shear  stress  distribution  in  longitudinal  direc- 

tion, upper  shelter,  30°  oblique  front-on  incidence.  Calculation  ID  3-50 

4- 1  Bending  moment  coefficients  at  center  of  closure,  front-on 

incidence,  comparison  of  finite  element  and  plate  theory  results  4-7 

4-2a  Velocity/time  history  for  a point  at  center  of  closure, 

longitudinal  component,  front-on  incidence.  Calculation  IB  4-8 

4-2b  Velocity/time  history  for  a point  at  center  of  closure, 

vertical  component,  front-on  incidence.  Calculation  IB  '-9 


8 


LIST  OF  FIGURES  (CONTINUED) 


Figure  No.  Page 

4-3a  Comparison  of  deformation  patterns  at  t = 10  msec.,  exaggerated, 

front-on  incidence  4-10 

4- 3b  Comparison  of  deformation  patterns  at  t = 25  msec. , exaggerated, 

front-on  incidence  4-11 

4-4a  Axial  stress  time  histories,  tubular  portion  behind  headworks, 

modified  Erode,  front-on  incidence  4-12 

4-4b  Axial  stress  time  histories,  tubular  portion  half  tube  diameter 

behind  headworks,  modified  Erode,  front-on  incidence  4-13 

4- 5  Comparison  of  breathing  mode  (n  = 0)  time  histories,  side-on 

incidence 4-14 

5- 1  S4  shelter  configuration 5-10 

5-2  S4  headworks 5-11 

5-3a  Three-dimensional  FE  model  of  S4  used  in  Phase  2 analysis, 

plane  of  symmetry  assumed 5-' 2 

5-3b  S4  soil  island  (front  view)  5-13 

5- 3c  Cross-section  of  S4  soil  island  and  shelter  perpendicular 

to  axis  of  the  tube 5-14 

5-3d  S4  headwork  model 5-15 

5-4a  Soil  model  in  uniaxial  strain,  axial  stress  versus  axial  strain  5-16 

5-4b  Soil  model  in  uniaxial  strain,  stress  difference  versus  pressure  5-17 

5-4c  Comparison  of  backfill  models,  uniaxial  stress/strain  curves  5-18 

5-5  Comparison  of  overpressure  waveforms  at  front  face,  front-on 

incidence.  Calculations  2A,  2E  and  2C 5-19 

5-6  Deformation  mode  of  headworks  at  123  msec,  10  msec  after  airbltst 

first  arrival  (113  msec),  head-on  incidence.  Calculation  2A  5-20 

5-7  Deformation  mode  of  backplate,  vertical  centerline.  Calculation  2A  5-21 

5-8  Deformation  mode  of  backplate,  horizontal  centerline.  Calculation  2A  . . . . 5-22 

5-9a  Development  of  inelastic  zones  in  backplate.  Calculation  2A. 

Arrival  time  of  airblast  is  113  msec 5-23 

5-9b  Development  of  Inelastic  zone  in  backplate.  Calculation  2A  (continued). 

Arrival  time  of  airblast  is  113  msec 5-24 

5-IOa  Development  of  inelastic  zone,  in  closure  concrete,  Calculation  2A. 

Arrival  time  of  airblast  is  113  msec 5-25 

5-lOb  Development  of  inelastic  zone,  in  closure  concrete,  Calculation  2A 

(continued).  Arrival  time  of  airblast  is  113  msec 5-26 

5-lOc  Development  of  inelastic  zone,  in  closure  concrete.  Calculation  2A 

(continued).  Arrival  time  of  airblast  is  113  msec 5-27 

5-lOd  Development  of  inelastic  zone,  in  closure  concrete.  Calculation  2A 

(continued).  Arrival  time  of  airblast  Is  113  msec 5-28 


9 


LIST  OF  FIGURES  (CONTINUED) 


Figure  No.  Page 

5-lOe  Development  of  inelastic  zone,  in  closure  concrete.  Calculation  2A 

(continued).  Arrival  time  of  airblast  is  113  msec 5-29 

5-11  Velocity/time  histories  of  closure  eg.  Calculation  2A  5-30 

5-12  Typical  stress  invariant  plot  for  points  in  closure  concrete  near 

support.  Calculation  2A  5-31 

5-13  Typical  stress  invariant  plot  for  points  in  center  of  closure 

concrete.  Calculation  2A  5-32 

5-14  Typical  stress  invariant  plot  for  points  in  backplate  near  support. 

Calculation  2A 5-33 

5-15  Typical  stress  invariant  plot  for  points  in  backplate  unsupported. 

Calculation  2A 5-34 

5-16a  Stress  invariant  plot  for  upper  frame  corner,  inner  steel  lining. 

Calculation  2A 5-35 

5-16b  Stress  invariant  plot  for  upper  frame  corner,  outer  steel  lining. 

Calculation  2A 5-36 

5-17a  Stress  Invariant  plot  for  lower  frame  corner.  Inner  steel  lining. 

Calculation  2A 5-37 

5-17b  Stress  Invariant  plot  for  lower  frame  corner,  outer  steel  lining. 

Calculation  2A 5-38 

5-18  Typical  stress  invariant  plot  for  concrete  in  upper  frame  corner. 

Calculation  2A 5-39 

5-19  Typical  stress  invariant  plot  for  concrete  In  lower  frame  corner. 

Calculation  2A 5-40 

5-20  Typical  stress  invariant  plot  for  concrete  In  frame  gussets. 

Calculation  2A 5-41 

5-21  Typical  stress  invariant  plot  for  points  In  headworks.  Calculation  2A  ...  5-42 

5-22  Typical  stress  Invariant  plot  for  points  In  the  tubular  section 

near  headworks,  crown  and  springline.  Calculation  2A  5-43 

5-23  Longitudinal  stress/time  histories  for  points  along  upper  shelter. 

Calculation  2A 5-44 

5-24  Longitudinal  stress/time  histories  for  points  along  lower  shelter, 

Calcul.T*’*"’!  2A 5-45 

5-25  Bearing  pressure  at  inner  support.  Calculation  2A  5-46 

5-26  Bearing  pressure  at  outer  support.  Calculation  2A  5-47 

5-27  Longitudinal  stress/time  histories,  tube  sections  behind  headworks. 

Calculation  2A 5-48 

5-28  Hoop  stress/time  histories  for  tube  section  behind  headworks. 

Calculation  2A 5-49 

5-29  Hoop  stress/time  histories  for  tube  section  half  tube  diameter 

behind  headworks.  Calculation  2A  5-50 


LIST  OF  FIGURES  (CONTINUED) 


1 

f 

\ 


A 

( 


( 


I 


FlRure  No.  Page 

5-30  Deformation  mode  of  cross-section  of  tube  (for  illustration  purpose 

only).  Section  A is  immediately  behind  headworks.  Section  B is  half- 

tube  diameter  behind  headworks  5-51 

5-31  Typical  strain  invariant  plot  for  points  in  backplate  near  support. 

Calculation  2A 5-52 

5-32  Typical  strain  invariant  plot  for  points  in  closure  concrete  near 

support.  Calculation  2A  5-53 

5-33  Typical  strain/time  histories  for  concrete  in  closure,  near  support. 

Calculation  2A 5-54 

5-34  Typical  straln/time  histories  for  concrete  in  closure,  near  center. 

Calculation  2A 5-55 

5-35  Deformation  modes  of  headworks,  head-on  incidence.  Calculation  2B  5-56 

5-36  Deformation  modes  of  backplate  vertical  centerline.  Calculation  2B  5-57 

5-37  Deformation  modes  of  backplate  horizontal  centerline.  Calculation  2B  . . . . 5-58 

5-38  Velocity/time  histories  of  closure  eg.  Calculation  2B  5-59 

5-39a  Typical  stress  Invariant  plot  for  points  at  closure  concrete,  near 

support.  Calculation  2B  5-60 

5- 39b  Typical  stress  invariant  plot  for  points  at  closure  concrete, 

center.  Calculation  2B  5-61 

5-40a  Typical  stress  invariant  plot  for  points  in  backplate,  near  support. 

Calculation  2B  ......  5-62 

5-40b  Typical  stress  invariant  plot  for  points  in  backplate,  near  support. 

Calculation  2B 5-63 

5-40c  Typical  stress  invariant  plot  for  points  in  backplate,  center. 

Calculation  2B 5-64 

5-41  Typical  stress  Invariant  plot  for  concrete  In  upper  frame 

corner.  Calculation  2B  5-65 

5-42  Typical  stress  invariant  plot  for  concrete  in  lower  frame 

corner.  Calculation  2B  5-66 

5-43  Stress  invariant  plot  for  steel  in  upper  frame  corner.  Calculation  2B  ...  5-67 

5-44  Stress  invariant  plot  for  steel  In  lower  frame  corner.  Calculation  2B  ...  5-68 

5-45  Bearing  load  time  histories.  Calculation  2B 5-69 

5-46  Longitudinal  stress/time  histories,  tube  sections  behind  headworks. 

Calculation  2B 5-/U 

5-4 7a  Hoop  stress/time  histories,  tube  section  immediately  behind  headworks. 

Calculation  2B 5-71 

5-47b.  Ihtop  stress/t Ime  histories,  tube  section  half-tube  diameter  behind 

headworks.  Calculation  2B  5-72 

5-48  Deformation  modes  of  headworks,  head-on  incidence.  Calculation  2C  5-73 

5-49  Deformation  modes  of  backplate  vertical  centerline.  Calculation  2C  5-74 


I 


i 

I 

I 


1 1 


!^s 


LIST  OF  FIGURES  (CONTINUED) 


Figure  No. 


Deformation  modes  of  backplate  horizontal  centerline.  Calculation  2C  ...  . 5-75 

Velocity /time  histories  of  closure  eg.  Calculation  2C  5-76 

Typical  stress  invariant-plot  for  points  in  closure  concrete,  near 

support.  Calculation  2C  5-77 

Typical  stress  invariant  plot  for  points  in  closure  concrete,  center. 

Calculation  2C 5-78 

Typical  stress  invariant  plot  for  points  in  backplate,  near  support. 

Calculation  2C 5-79 

Typical  stress  invariant  plot  for  points  in  backplate,  near  center. 

Calculation  2C 5-80 

Typical  stress  invariant  plot  for  concrete  in  upper  frame  corner. 

Calculation  2C 5-81 

Typical  stress  invariant  plot  for  concrete  in  lower  frame  corner. 

Calculation  2C 5-82 

Typical  stress  invariant  plot  for  steel  in  upper  frame  corner. 

Calculation  2C 5-83 

Typical  stress  invariant  plot  for  steel  in  lower  frame  corner, 

Calculation  2C 5-84 

Typical  stress  invariant  plot  for  concrete  in  tubular  section, 

crown  and  springline.  Calculation  2C 5-85 

Typical  stress  invariant  plot  for  concrete  in  tubular  section, 

invert.  Calculation  2C  5-86 

Bearing  load/time  histories.  Calculation  2C  5-87 

Axial  stress/time  histories  for  points  in  upper  shelter.  Calculation  2C  . . 5-88 
Axial  stress/time  histories  for  points  in  lower  shelter.  Calculation  2C  . . 5-89 

Comparison  of  velocity/time  histories  of  closure  eg  6-7 

Model  for  TRANAL  check  problem A-6 

Comparison  of  membrane  stress  in  steel  liner  A-7 

Comparison  of  radial  displacements  in  steel  liner  A-8 


LIST  OF  TABLES 


Table  No.  Page 

3-1  Shelter  "Ring"  In  Vacuo  Frequencies  (Hz)  3-51 

3-2  Elastic  Properties  of  Soil  and  Structure 3-52 

5-1  Material  Parameters  for  S4  Model 5-90 


13 


SECTION  I 


INTRODUCTION 

The  shelter  concept  Is  one  of  the  many  multiple  aim  point  (MAP)  basing  mode  con- 
cepts under  consideration  for  deployment  of  the  M-X  advanced  ballistic  missile  system.  In 
the  shelter  based  mode,  missiles  are  moved  by  transporter/launcher  vehicles  among  widely 
dispersed  shelters  as  compared  to  the  current  Minuteman  ballistic  missile  force  which  is 
based  in  hardened  silos  arranged  in  clusters  or  wings.  Random  shifting  of  the  M-X  missiles 
in  the  MAP  shelters  would  Increase  their  survivability  by  complicating  the  targeting  options 
of  an  enemy  attack. 

Although  many  different  shelter  configurations  are  under  consideration  at  this 
time,  there  are  some  structural  characteristics  which  are  common  to  all  configurations  and 
a structure  is  said  to  be  shelter-like  if  it  displays  these  characteristics:  if  it  is  buried, 
it  is  burled  with  the  tube  horizontal  and  the  front  face  vertical  or  near  vertical.  The 
shelter  may  house  the  transporter /carrier  in  addition  to  the  missile  and  in  that  case  some 
portion  of  it  must  be  large  enough  to  accommodate  the  vehicle.  This  larger  portion  is 
joined  to  the  tube  portion  through  a transition  section.  Access  to  the  shelter  is  through 
one  end  (the  front  end)  which  is  normally  guarded  by  a door  or  closure,  while  the  other 
end  is  closed.  To  enable  the  transporter/carrier  to  move  in  and  out  of  the  shelter  with 
some  speed  and  efficiency,  a key  to  the  shelter  basing  concept,  a driveway  leading  to  the 
door  and  front  end  is  necessary.  The  driveway  is  below  grade  if  the  shelter  is  below 
grade. 

The  large  front-end  portion  of  the  shelter  (called  the  headworks)  is  usually 
rectangular  in  shape  and  has  a rectangular  cavity  to  house  the  carrier/vehicle.  Hence, 
the  closure  is  also  rectangular.  When  the  missile  alone  is  to  be  accommodated,  the  rec- 
tangular front  end  is  no  longer  a requirement.  However,  some  form  of  a front  end, 
more  massive  than  the  shelter  tube  proper  is  still  necessary  to  house  the  closure 
envisioned.  In  this  case,  the  front  end  (and  the  closure)  can  also  be  circular  in 
order  to  mesh  in  more  naturally  with  the  tubular  section.  The  transition  between  front 
end  and  tube  may  be  less  abrupt  as  in  the  rectangular  headworks  case.  This  study  deals 
with  the  more  complex  analysis  problem  of  the  rectangular  headworks  and  closure  but  the 
results  may  also  be  applied  to  a circular  headworks  and  closure  configuration. 

Whereas  the  shelter  is  protected  on  its  top,  flanks  and  aft-end  from  the  hostile 
environment  by  the  soil  medium  in  which  it  is  burled,  its  front  is  necessarily  exposed  to 
nuclear  alrblasts  because  of  the  open  accessway.  What  compounds  this  situation  is  the  fact 
that  the  front  end  la  vertical  or  near  vertical.  When  the  alrblast  impinges  on  it  head-on. 
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the  overpressure  is  greatly  magnif led  due  to  the  reflection  phenomenon.  A central  prob- 
lem In  the  design  of  a shelter-llke  structure  is  the  design  of  the  closure  and  frame 
(front  face)  which  must  withstand  this  tremendous  load.  Dynamic  analyses  of  the  structure 
are  dominated  by  the  analysis  of  the  front-end  response  and  its  effects  on  the  remainder 
of  the  structure. 

The  present  study  Is  conducted  mainly  to  seek  answers  to  the  following  questions: 

• What  are  the  deformation  modes  of  the  shelter  under  different  loading  condi- 
tions (slde-on,  head-on,  oblique  front-on),  and  how  sensitive  are  these  deformation  modes 
to  variations  in  the  overpressure  wave  forms? 

• What  are  the  important  features  of  the  shelter  which  must  be  Included  in  the 
analysis  model,  and  how  should  these  features  be  modeled? 

• What  Is  the  correct  use  of  an  elastic  analysis  of  the  shelter  structure?  What 
pertinent  information  can  be  and  cannot  be  obtained  from  such  an  elastic  analysis? 

• What  are  the  main  effects  of  inelasticity  and  how  important  is  an  inelastic 
analysis  in  the  analysis  of  the  shelter? 

• What  is  the  role  of  ground  motion  in  the  shelter  analysis?  Is  the  soil  island 
technique  necessary  for  a complete  analysis  of  the  shelter  structural  response? 

Through  answers  to  these  questions  certain  guidelines  emerge,  culminating  in  a 
methodology  that  can  be  followed  in  performing  future  dynamic  analyses  of  a shelter-llke 
structure  in  a nuclear  attack  environment.  This  methodology  constitutes  the  main  result 
of  the  present  study. 

A byproduct  of  this  effort  is  the  wealth  of  quantitative  results  computed  for 
certain  prototype  shelters  which  have  served  as  case  studies.  The  data  may  be  used  to 
aid  in  the  evaluation  and  validation  of  these  shelter  configurations  (SI  and  S4).  Whereas 
numerical  data  tend  to  become  obsolete  as  the  design  is  changed,  it  is  believed  that  the 
methodology  applicable  to  a class  of  shelter-llke  structures  may  have  longer  lasting  value 
to  the  study  of  protective  structures. 
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SECTION  II 
SCOPE 

The  scope  of  this  effort  as  originally  outlined  in  Reference  1 has  undergone 
several  modifications  during  the  course  of  the  study  in  response  to  different  project  re- 
quirements. Due  to  the  lack  of  overpressure  data  and  the  need  to  support  SAMSO's  closure 
design  validation  effort,  the  study  has  deviated  from  a research  type  to  include  limited 
applications;  in  particular  it  includes  some  detailed  analyses  of  the  S4  test  configuration. 

The  scope  presented  herein  reflects  these  changes,  and  can  be  divided  into  three 
phases:  Phase  1 which  identifies  the  elastic  structural  deformation  modes  and  their  sen- 
sitivity to  variations  in  the  applied  airblast  wave  forms.  Phase  2 which  provides  response 
definitions  for  a particular  shelter  configuration  subjected  to  a particular  airblast  load 
in  order  to  show  the  effect  of  Inelasticity  and  modeling  techniques,  and  Phase  3 which 
comprises  the  work  performed  in  support  of  the  S4  test  and  closure  design  evaluation. 

Phases  1 and  2 will  be  described  herein,  and  because  of  its  somewhat  different  purpose 
Phase  3 will  be  described  in  a subsequent  report  (Reference  2). 

The  shelter-like  structure  is  assumed  to  be  located  at  a generic  site,  at  the  ^ 

600  psi  range  from  a 1 MT  nuclear  explosion,  or  about  1,850  feet  from  ground  zero.  The 

airblast  wave  front  can  impinge  the  shelter/berm  configuration  head-on,  at  an  oblique  ^ 

angle,  side-on  or  rear-on.  This  study  addresses  only  the  head-on,  oblique  (30°)  and  side- 
on  cases  with  emphasis  on  head-on  incidence. 

2.1.  PHASE  ONE 

As  a first  step,  elastic  calculations  were  performed  to  identify  basic  struc- 
tural response  modes  of  the  MAP  shelter.  These  calculations  are  both  two-dimensional,  in 
which  the  cylindrical  portions  of  the  shelter  are  modeled  as  a ring,  and  three-dimensional 
in  which  the  shelter  is  modeled  in  its  entirety.  In  the  head-on  configuration,  where  the  ; 

airblast  impinges  on  the  front  face  head-on,  there  is  a vertical  symmetry  plane  containing 

the  longitudinal  axis;  consequently,  only  half  of  the  shelter/berm  ensemble  needs  to  be  ■ 

modeled.  When  the  airblast  approaches  the  front  face  at  an  oblique  angle  (other  than  head-  i 

on),  this  symmetry  is  lost  and  the  entire  shelter/berm  ensemble  is  modeled.  All  shelter  * 

I 

structures  studied  in  Phase  I correspond  to  the  HAVE  HOST  SI  configuration  (Reference  3).  | 

I 

The  calculations  were  performed  using  GENSAP,  an  implicit  elastic  three- 
dimensional  finite  element  code.  The  major  shelter/loading  configurations  studied  are 
as  follows: 
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Calculation  lA  has  a slde-on  loading  of  the  cylindrical  portion  of  the  shelter 
from  an  air  overpressure  which  passes  over  the  berm,  with  subsequent  modification  in  both 
pressure  and  time  characteristics  due  to  reflections  and  rarefactions  by  the  berm.  This 
calculation  is  two  dimensional. 

Calculation  IB  has  a front-on  loading  in  which  the  air  overpressure  directly 
loads  the  front-face  of  the  shelter  and  then  progresses  over  the  berm  in  a direction 
parallel  to  the  longitudinal  axis  of  the  structure.  Because  of  the  symmetry  in  structure 
and  overpressure  with  respect  to  the  center  plane,  only  half  of  the  physical  configuration 
is  included  in  the  three-dimensional  calculation. 

Calculation  1C  is  similar  to  Calculation  IB  with  one  exception.  The  monolithic 
front  face  used  in  IB  is  modified  to  reflect  the  existence  of  a gap  between  the  closure 
and  frame.  Air  overpressure  which  acts  on  the  front  face  is  allowed  to  penetrate  into 
this  gap.  By  comparing  the  results  of  Calculations  IB  and  1C,  the  effect  of  the  gap  on 
structural  response  is  identified. 

Calculation  ID  has  a front-on  loading  in  which  the  direction  of  airblast  propaga- 
tion is  inclined  at  30°  to  the  normal  to  the  front  face  which  is  modeled  as  a monolithic 
structure  (i.e.,  no  gap).  By  comparing  the  results  of  Calculations  IB  and  ID,  the  effect 
of  oblique  incidence  is  ascertained. 

In  addition  to  identifying  the  basic  response  modes  to  loading  conditions,  it 
was  the  Intention  of  the  study  to  investigate  the  sensitivity  of  these  response  modes  to 
variations  in  the  overpressure  wave  forms.  In  particular,  the  following  three  types  of  wave 
forms  were  considered  for  use  in  the  study: 

a.  Free-air  wave  form  from  a 1 MT  surface  burst. 

b.  Precursed  wave  form  from  a 1 MT  surface  burst  with  thermal  layer. 

c.  Precursed  wave  form  from  a 1 MT  surface  burst  in  dusty  air. 

However,  since  airblast  data  available  from  the  Air  Force  Weapons  Laboratory  (AFWL) 
were  limited,  all  Phase  1 calculations  used  basically  the  free  air  wave  form  from  a 1 MT 
surface  burst  (References  4 and  5).  Consequently,  the  question  of  how  sensitive  the  struc- 
tural response  is  to  variations  in  overpressure  wave  forms  remains  unanswered.  To  gain  some 
Insight  into  this  sensitivity,  however,  results  of  Calculations  lA  to  1C  have  been  compared 
with  results  obtained  in  some  checkout  calculations,  which  use  the  Erode  overpressure  wave 
form  for  flat  ground,  modified  in  some  cases  to  account  for  reflections  at  the  shelter  front 
face. 
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All  Phase  1 calculations  are  made  using  alrblast  loadings  only  and  only  Che 
short-term  response  (e.g.,  somewhat  greater  times  than  two  transit  times  of  the  shock  wave 
across  Che  cylindrical  section  of  the  shelter  In  the  two-dimensional,  slde-on  case  and  i 

somewhat  greater  times  Chan  the  times  of  peak  stress  In  the  headworks  and  transition  sec- 
tion In  Che  three-dimensional  cases)  Is  obtained. 

2.2.  PHASE  TWO 

In  this  phase  of  the  study,  Che  Inelastic  structural  response  Is  analyzed.  The 
behavior  of  the  structure  under  alrblast  loading  Identified  In  the  first  phase  Is  quantified,  | 

and  deformation  modes  due  to  Inelastic  material  behavior  are  defined.  These  Inelastic  cal-  | 

culatlons  are  three-dimensional.  Where  symmetry  exists,  as  In  the  head-on  loading  case, 
only  half  of  the  shelter/berm  ensemble  Is  modeled.  All  shelter  structures  studied  In  this 
phase  correspond  to  the  HAVE  HOST  S4  configuration,  or  Its  replacement  the  HEST  S4  con- 
figuration (Reference  6).  Phase  2 study  consists  of  the  following  three-dimensional  cal- 
culations: 

Calculation  2A  has  a front-on  loading  condition.  The  concrete  In  the  structure 
is  modeled  as  an  elastic  perfectly  plastic  material  and  the  back  and  side  plates  in  the 
closure  are  assigned  a yield  strength  corresponding  to  that  of  A36  steel.  The  steel  plates 
are  modeled  by  two  elements  across  the  plate  thickness.  The  backfill  model  corresponds  to 
the  medium  backfill  of  HAVE  HOST  (Reference  7).  | 

Calculation  2B  Is  similar  to  Calculation  2A  except  that  the  concrete  is  the  only 
Inelastic  material  assumed.  The  elastic  back  and  side  plates  are  modeled  by  only  one  ele- 
ment across  its  thickness.  The  backfill  has  the  elastic  unloading  properties. 

Calculation  2C  Is  similar  to  Calculation  2B,  but  with  all  materials  (concrete, 
steel,  backfill,  soil)  elastic. 

Whereas  the  shelter  structure  analyzed  in  Phase  1 corresponds  to  the  SI  configura-  | 

tion,  the  structure  analyzed  in  Phase  2 corresponds  to  the  S4  configuration.  This  change  1 

, In  configuration  is  made  to  support  the  closure  design  validation  program;  It  tends  to  make  1 

{ correlation  of  the  results  obtained  In  Phases  1 and  2 more  difficult.  ^ 

: The  air  overpressure  used  as  Input  to  the  Phase  2 calculations  Is  provided  by 

[ the  AFWL  (Reference  8),  based  partly  on  test  data,  partly  on  design  handbooks  and  partly 

' on  two-dimensional  numerical  calculations. 
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All  Phase  2 calculations  are  made  using  Weidlinger  Associates'  TRANAL  code,  an 
explicit  inelastic  three-dimensional  dynamic  finite  element  program  (Reference  10).  To 
accommodate  the  rather  complicated  geometry  and  applied  load  distribution  of  the  S4  con- 
figuration, extensive  modifications  are  made  to  the  code  and  numerous  subroutines  written 
with  most  of  the  changes  occurring  in  the  input  module,  element  table  processing,  coordinate 
and  applied  load  processing.  A brief  summary  of  these  changes  and  a description  of  the  add- 
ed code  capability  are  given  in  Appendix  A. 

Results  of  Phases  1 and  2 calculations  are  described  in  Sections  III  and  V, 
respectively.  Although  it  is  not  the  purpose  of  this  study  to  examine  structural  response 
of  Che  shelter  in  detail,  some  description  of  the  structural  response  is  necessary  in  order 
to  lead  into  the  main  discussion  of  the  study  results.  Accordingly,  the  description  in 
Sections  III  and  V will  be  brief  and  is  presented  on  a calculation  basis,  as  if  each  cal- 
culation is  performed  alone. 

These  results  are  compared  and  contrasted  in  Section  IV  for  Phase  1 and  in  Sec- 
tion VI  for  Phase  2.  Through  such  comparisons  it  is  possible  to  (a)  identify  the  shelter 
deformation  modes  and  their  relationship  to  loading  conditions  and  airblast  characteristics; 
(b)  isolate  vulnerable  shelter  elements  and  their  modeling  requirements;  and  (c)  clarify 
the  effect  of  inelasticity  on  shelter  response  and  hence  the  importance  of  an  inelastic 
analysis.  The  conclusions  presented  in  Section  VII  lead  naturally  to  a procedure  for  the 
dynamic  analysis  of  the  shelter  structure,  the  major  goal  of  the  present  study. 
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PHASE  1 RESULTS 

Results  of  calculations  performed  in  Phase  1 will  be  described  in  this  section. 

Data  presented  consist  of  time  histories  of  sti  ess/strain  in  the  elements  and  velocity/ 
displacement  of  the  nodes.  A brief  description  of  the  finite  element  model,  the  air  over- 
pressure and  ground  motion  input  will  be  given  where  appropriate. 

All  Phase  1 calculations  are  based  on  a prototype  shelter  corresponding  to  the 
HAVE  HOST  SI  test  structure  as  shown  In  Figure  3-1  (Reference  3).  The  headworks  is  sym- 
metric on  two  planes  and  the  tube  axis  is  coincident  with  the  headworks  axis.  The  closure, 
a thick  concrete  slab,  is  simply  bolted  to  the  bearing  area.  In  this  configuration  the 
front  face  and  part  of  the  sides  and  top  portion  of  the  headworks  are  exposed  to  direct  air- 
blast  load.  There  are  no  headwalls. 

3.1.  CALCULATION  lA— TITO-DIMENSIONAL  ELASTIC  CALCULATION  WITH  SIDE-ON  INCIDENCE 

The  two-dimensional  calculation  considers  deformation  in  a plane  perpendicular 
to  the  shelter  longitudinal  axis  and  the  shelter  is  represented  by  a circular  cross  sec- 
tion of  the  cylindrical  portion,  l.e.  , a ring.  No  displacement  is  allowed  in  the  direc- 
tion of  the  shelter  longitudinal  axis. 

Model  and  Airblast  Load — the  plane  strain  model  is  shown  in  Figure  3-2.  It 
consists  of  plane  strain  quadrilateral  finite  elements  (approximately  1,700  in  soil  and 
70  in  the  structure).  The  refinement  in  the  structure  is  determined  from  a preliminary 
finite  element  calculation  of  natural  mode  shapes  and  frequencies  of  the  structure  in  vacuo 
which  are  compared  with  its  analytical  solutions,  such  as  shown  in  Table  3-1.  Improved 
agreement  for  higher  modes  could  be  obtained  by  refining  the  grid  further,  but  this  is  deem- 
ed unnecessary  for  the  purpose  of  the  present  study. 

The  elastic  material  properties  for  the  native  soil,  berm  and  backfill  are  given 
in  Table  3-2.  The  grid  sizes  of  the  individual  finite  elements  are  selected  so  as  to  pro- 
vide undistorted  components  of  airblast-induced  ground  motion  up  to  frequencies  of  40  Hz 
in  the  native  soil  and  20  to  25  Hz  for  the  berm  and  backfill.  The  overall  size  of  the  model 
is  selected  to  give  approximately  80  msec,  of  useful  structural  response  time,  or  about  two 
cycles  of  the  ovaling  (n  = 2)  mode. 

The  airblast  used  in  this  calculation  is  generated  by  the  AFWL  HULL  code,  designat- 
ed HULL  19.8042  (Reference  5).  It  is  generated  for  free-air,  but  takes  into  account  reflec- 
tions due  to  the  existence  of  the  berm.  Later  on  in  the  discussion  reference  will  be  made 
to  the  Erode  description  which  is  the  classical  nuclear  wave  form  for  a flat  rigid  ground. 


Because  of  the  limited  data  on  air  overpressure  which  are  available  a direct  investiga- 
tion of  the  effect  of  variations  in  overpressure  wave  form  is  not  feasible.  By  comparing 
response  results  obtained  using  the  HULL  19.8042  and  the  Brode  wave  forms,  even  though  the 
latter  is  not  compatible  with  the  berm,  some  insight  into  the  wave  form  effect  can  be 
obtained. 

A comparison  of  the  peak  overpressure  profile  over  the  berm  surface  for  the  two 
alrblast  descriptions  is  given  in  Figure  3-3.  To  obtain  an  understanding  of  the  differ- 
ences in  temporal  characteristics  between  the  two  alrblast  descriptions,  the  overpressure 
time  histories  at  selected  points  on  the  berm  surface  are  scaled  to  a constant  peak  of  unit 
magnitude  and  the  scaled  histories  compared  in  Figure  3-4.  Note  tK'at  the  temporal  charac- 
teristics of  the  two  alrblast  wave  forms  are  quite  similar  for  all  points  on  the  berm  with 
possibly  the  exception  of  the  side  away  from  GZ  which  is  considered  to  have  very  minor 
effect  on  structural  response.  Hence,  it  is  fair  to  say  that  any  sensitivity  of  structural 
response  to  wave  form  is  probably  induced  by  differences  in  peak  overpressures  as  illus- 
trated in  Figure  3-3,  the  more  significant  of  which  occurs  along  the  blast  side  and  also 
at  the  top  of  the  berm. 

Deformation  Mode — the  deformation  pattern  of  the  structure  is  mainly  that  of  the 
breathing  and  ovallng  (n  = 2)  mode.  This  is  illustrated  in  Figure  3-5  which  shows  the 
radial  displacement  (not  to  same  scale  as  the  ring  radius)  as  a function  of  the  circumfer- 
ential distance  along  the  ring  at  different  action  times.  The  dotted  line  corresponds  to 
the  breathing  mode  contribution.  It  can  be  seen  that  aside  from  the  breathing  mode  the 
deformation  is  dominated  by  the  n = 2 mode;  the  finding  is  also  confirmed  by  a Fourier 
analysis  of  the  deformation  pattern,  shown  in  Figure  3-6. 

Stresses — the  stress/time  histories  for  elements  in  the  structure  are  consistent 
with  the  deformation  pattern  just  described.  The  circumferential  stresses  for  the  three 
elements  which  make  up  the  thickness  of  the  ring  can  be  used  to  estimate  the  amount  of 
bending  and  hoop  compression  which  exist  in  the  structure.  Figure  3-7  shows  the  circum- 
ferential stresses  at  the  crown  (12  o'clock  position)  for  the  Innermost  and  outermost  ele- 
ments; these  can  be  averaged  to  yield  the  hoop  compression  assuming  plane  sections  remain 
plane.  The  difference  between  the  circumferential  stresses  and  hoop  compression  is  propor- 
tional to  the  bending  moment.  Similar  time  histories  are  obtained  for  points  at  the  bottom 
of  the  structure  (6  o'clock  position)  and  at  the  springlines  (9  o'clock  position  and 
3 o'clock  position). 

The  maximum  bending  moment  thus  obtained  is  approximately  2.4  x 10*  in. -lb. /in., 
and  the  maximum  membrane  stresses  due  to  bending,  i.e.,  bending  stress  at  the  extreme 
fibre,  is  11,000  psl. 


Soil-Structure  Interaction — the  first  wave  to  reach  the  structure  is  from  the  top 
of  the  berm  which  is  about  10  feet  of  overburden.  The  structure  then  senses  the  airblast 

load  coming  from  the  windward  side  of  the  berm  and  finally  that  from  the  leeward  side  of 

the  berm.  There  is  some  evidence  of  outrunning  ground  motion  in  the  native  soil  underly- 
ing the  backfill  but  its  effect  does  not  appear  significant. 

The  initial  airblast  Induced  ground  motion  wave  front  is  reflected  at  the  crown 
of  the  shelter  and  again  reflected  at  the  free-air  boundary.  Consequently,  significant 
ringing  (vibration)  is  observed  in  the  soil  motion  time  histories  (as  Figure  3-8  shows). 
This  effect  which  is  prominent  in  an  elastic  calculation  is  expected  to  diminish  greatly 
when  more  realistic  nonlinear  soil  properties  are  used  to  model  hysteresis  and  other  dis- 
sipative mechanisms  inherent  in  the  backfill. 

Another  soil-structure  Interaction  phenomenon  manifests  itself  at  depths  lower 
than  the  crown,  where  the  first  loading  of  the  soil  is  Induced  by  waves  propagated  through 
the  shelter  and  out  into  the  soil  (an  outrunning  effect).  The  structure  is  flattened  by 
the  ground  shock  at  the  crown  and  pushes  out  at  the  soil  at  the  springline.  This  effect 
can  be  observed  before  the  arrival  of  the  direct  airblast-induced  ground  motion.  Figure 
3-9  shows  the  stress-time  histories  for  a soil  element  left  on  the  springline;  the  hori- 
zontal stress  (O^)  Induced  by  motion  of  the  structure  precedes  the  main  airblast  load 
which  is  almost  all  vertical  (a  ).  Similarly,  Figure  3-10  Illustrates  this  phenomenon 
for  a soil  element  to  the  right  of  the  springline.  Note  the  synchronization  of  the 
pulse  with  that  of  Figure  3-9  and  the  main  pulse  from  the  overpressure  load  acting  at  the 
top  and  sides  of  the  beam  arriving  later  in  time. 

3.2.  CALCULATION  IB— THREE-DIMENSIONAL  ELASTIC  CALCULATION  WITH  HEAD-ON  INCIDENCE 

This  three-dimensional  calculation  considers  deformation  symmetric  to  the  center 
plane  of  symmetry  of  the  structure/berm  configuration.  Hence,  only  one  half  of  the  con- 
figuration needs  to  be  modeled.  The  alrblast-load  impinges  on  the  front  face,  progresses 
over  the  top  and  then  travels  down  the  berm  in  the  direction  of  the  longitudinal  axis  of 
the  shelter. 

Model  and  Airblast  Load — a view  of  the  three-dimensional  finite  element  model 
which  consists  of  eight  node  (brick)  elements  is  shown  in  Figure  3-11.  The  model  consists 
of  approximately  1,000  elements  in  soil  and  600  elements  in  the  structure  and  covers  the 
full  length  of  the  structure  with  approximately  40  feet  of  soil  in  front  of,  20  feet  below 
and  to  the  side  of  the  structure.  A cross-sectional  view  of  the  model  is  shown  in  Figure 
3-12  and  Figure  3-13  gives  another  view  of  the  structure,  with  soil  elements  removed.  Note 
that  special  attention  has  been  given  to  the  transition  region,  where  the  main  headworks 
change  from  a rectangular  box  to  a circular  cylinder. 
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In  this  calculation  the  front  face  is  modeled  simply  as  a monolithic  plate,  i.e., 
no  distinction  is  made  between  the  closure  and  its  frame,  and  the  gap  between  them  is 

ignored.  The  gap  though  small,  allows  the  high  reflected  overpressure  which  exists  on  the 

front  face  a passage  inside  the  frame  and  is  shown  in  Section  3.3  to  have  a significant 
effect  on  the  behavior  of  the  frame  when  included  in  the  model.  Henceforth,  the  term 
no-gap  refers  to  the  situation  where  the  front  face  is  modeled  as  monolithic  and  gap  de- 
notes the  presence  of  the  gap  in  the  model. 

The  elastic  material  properties  for  the  native  soil,  berm,  backfill  and  struc- 
ture are  the  same  as  those  used  in  Calculation  lA.  the  elastic  side-on  two-dimensional 

calculation.  The  grid  sizes  on  the  individual  finite  elements  are  not  as  refined  as  would 

be  desirable.  For  instance,  the  cylindrical  portion  of  the  structure  is  modeled  by  two 
elements  across  the  thickness  of  the  ring  in  contrast  to  three  in  Calculation  lA  and  some 
of  the  soil  elements  below  and  to  the  side  of  the  structure  are  rather  large.  However, 
the  element  sizes  can  still  provide  undistorted  components  of  motion  up  to  frequencies  of 
160  Hz  in  the  structure  and  14  Hz  in  the  overburden. 

The  AFWL  HULL  19.8041  airblast  (Reference  4)  for  free  air  is  used  in  this  calcula- 
tion. Some  model  checkout  results  obtained  using  a modified  Erode  airblast  for  flat,  rigid 
ground  will  be  referred  to  in  a later  section.  This  corresponds  to  using  the  Erode  descrip- 
tion unchanged  everywhere  in  the  configuration  with  the  exception  of  the  front  face  where 
the  nominal  peak  overpressure  of  600  psi  is  scaled  to  5,400  psi,  the  peak  value  in  the  AFWL 
data.  A comparison  of  the  peak  overpressure  profiles  over  the  shelter/berm  surface  for  the 
two  airblast  descriptions  is  given  in  Figure  3-14.  Scaled  time  histories  at  selected  points 
on  the  structure/berm  surface  are  compared  in  Figure  3-15.  Differences  in  wave  form  are 
apparent,  especially  for  points  in  front  of  the  shelter  which  encounter  the  incoming  air- 
blast  and  then  the  reflected  wave  from  the  front  face.  A more  important  difference,  how- 
ever, is  in  the  impulse  or  the  area  under  the  curves.  Referring  to  Figure  3-15a  for  a point 
on  the  front  face,  the  Erode  curve  has  about  two  times  the  Impulse  as  the  HULL  19.8041  in 

the  first  20  msec,  alone.  Its  significance  will  be  apparent  when  the  response  results  are 

described. 

Both  airblast  descriptions  are  obtained  assuming  cylindrical  symmetry,  i.e., 
the  axis  of  symmetry  being  a line  at  the  burst  point  (ground  zero)  and  perpendicular  to 
the  ground.  Hence,  modifications  would  be  necessary  for  use  on  a three-dimensional  struc- 
ture. In  the  HULL  19.8041  description  the  side  of  the  front  face  is  assumed  to  have  the 

same  airblast  load  distribution  as  the  top  of  the  front  face.  In  the  modified  Erode  des- 

cription, all  points  on  the  front  face  are  assumed  loaded  uniformly.  When  the  gap  is 
modeled  (Section  3.3),  the  same  pressure  i.s  assumed  to  apply  within  the  gap.  Loads  at  the 
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exposed  side  and  top  surfaces  of  the  headworks  are  obtained  using  magnification  factors 
measured  in  the  HAVE  HOST  SI  test  (Reference  11). 

A comparison  of  the  peak  overpressure  distributions  on  the  front  face  for  these 
two  descriptions  is  given  in  Figure  3-16.  A portion  of  the  front  face  is  shielded  by  soil 
and  hence  is  not  subjected  to  direct  airblast  loading.  On  the  portion  which  is  exposed 
is  plotted  the  peak  overpressure  distributions.  The  dotted  line  represents  the  HULL 
19.8041  pressure  and  the  dashed  line  represents  the  modified  Erode  pressure  which  of 
course  is  constant  over  the  front  face.  The  volume  under  each  spatial  surface  is  the 
total  axial  load  acting  on  the  front  face  in  each  case.  The  resultant  force  acts  at  a 
point  in  the  front  face  which  is  offset  from  the  geometric  center  of  the  face.  The  pro- 
duct of  the  resultant  force  and  the  offset  gives  a measure  of  the  bending  moment  acting 
on  the  headwork  in  each  case. 

Deformation  Mode — the  deformation  pattern  of  the  loaded  structure  is  similar  to 
that  of  a beam-column  on  an  elastic  foundation.  The  tremendous  airblast  load  acting  di- 
rectly at  the  front  face  dictates  the  deformation  pattern  that  the  structure  follows.  The 
shelter  is  compressed  as  a column  and  because  of  the  load/resistance  imbalance  acting  at 
the  front,  l.e.,  offset  between  resultant  force  and  geometric  (stiffness)  centroid,  a cer- 
tain amount  of  bending  (rotation)  of  the  headworks  in  the  plane  of  symmetry  occurs;  the 
extent  of  this  bending  deformation  depends  on  the  magnitude  of  the  resultant  axial  load  and 
the  degree  of  load  imbalance  or  offset.  The  airblast  overpressure  acting  on  the  overburden 
tends  to  push  the  shelter  down  and  in  so  doing  contributes  somewhat  to  the  aforementioned 
longitudinal  bending.  These  two  actions,  however,  are  not  synchronous  since  one  is  effected 
by  direct  airblast  load  whereas  the  other  indirectly  by  induced  ground  motion  at  the  over- 
burden. The  deformation  mode  is  Illustrated  in  Figure  3-17  where  the  displacements  are 
exaggerated  for  illustration  purpose. 

Some  breathing  and  ovallng  deformation  in  planes  perpendicular  to  the  longitudinal 
axis  are  observed,  although  at  this  early  time,  they  have  not  been  fully  developed  compared 
to  the  axial  load  induced  deformations  (both  axial  compression  and  longitudinal  bending). 

Stresses — the  direct  load  on  the  front  face  produces  high  axial  stresses  in  the 
front  portion  of  the  structure  which  become  even  higher  as  the  stress  wave  propagates  into 
the  transition  region  and  finally  into  the  tubular  region  (see  Figures  3-18  and  3-19).  This 
follows  immediately  from  the  fact  that  the  load-bearing  cross-sectional  area  of  the  shelter 
decreases  as  one  moves  along  the  structure. 

For  a particular  cross-section  of  the  tube  the  bending  moments  and  stresses  in 
that  plane  can  be  estimated  in  exactly  the  same  manner  as  in  Calculation  lA. 
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Similarly,  the  bending  moment  responsible  for  longitudinal  bending  can  be 
estimated  using  axial  stress  time  histories  (Figures  3-20  and  3-21)  at  the  top  and 
bottom  of  the  shelter. 

Maximum  average  axial  compression  is  found  to  be  7,300  psl.  The  longitudinal 
bending  moment  is  estimated  to  be  3.4  x 10®  in. -lb.,  and  the  corresponding  extreme  fibre 
bending  stress  is  2,800  psl.  The  maximum  bending  and  compressive  stresses  occur  at  dif- 
ferent times. 

The  peak  pressure  exerted  by  the  closure  on  the  bearing  surface  has  a distribu- 
tion as  shown  in  Figure  3-22.  The  pressure  is  highest  at  the  top  bearing  ring  and  is  in 
general  higher  at  midspan.  The  average  peak  pressure  of  9,000  psl  is  slightly  less  than 
twice  the  peak  applied  load  whereas  the  ratio  of  the  closure  loaded  area  to  bearing  area 
is  2.5.  Since  the  front  face  of  the  headworks  is  monolithic  in  this  calculation  not  all 
of  the  closure  load  is  supported  by  the  bearing  area.  Some  support  is  afforded  by  the  sur- 
rounding frame  in  this  model. 

Soil-Structure  Interaction — at  early  times  the  berm  overburden  effectively  iso- 
lates the  structure  (tubular  section)  from  most  of  the  top  airblast  overpressure.  For 
instance  it  takes  a little  more  than  10  msec,  for  an  airblast  induced  ground  motion  to 
transverse  the  overburden  and  reach  the  crown  of  the  tubular  section  of  the  shelter. 

The  soil  shielding  the  lower  portion  of  the  shelter  front  causes  the  center  of 
pressure  to  become  offset  from  the  center  of  stiffness.  This  offset,  multiplied  by  the 
applied  load  generates  a bending  moment  and  is  the  imbalance  referred  to  in  earlier  para- 
graphs. 

3.3.  CALCULATION  1C— THREE-DIMENSIONAL  ELASTIC  CALCULATION  WITH  HEAD-ON  INCIDENCE  AND 
FRONT- FACE  GAP 

The  gap  refers  to  the  opening  between  sides  of  the  closure  and  the  headworks 
frame  surrounding  the  closure.  Despite  its  diminutive  dimension,  if  it  is  assumed  that 
the  high  reflection  pressure  which  acts  on  the  front  face  penetrates  into  the  gap  cavity 
as  well,  its  effect  on  stress  distribution  in  the  headworks  becomes  significant  as  shown 
by  results  of  this  calculation. 

Model  and  Airblast  Load — the  finite  element  model  used  in  Calcuation  IB  is  modi- 
fied to  Include  the  presence  of  a gap  between  closure  and  headworks.  This  is  achieved  by 
eliminating  the  connectivity  between  elements  on  the  two  sides  of  the  closure/frame  inter- 
face. All  other  aspects  of  the  model  remained  basically  unchanged. 
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The  alrblast  description  used  in  Calculation  IB,  namely  AFWL  HULL  19.8041,  is 
again  employed  in  order  to  assess  the  effect  of  the  gap  alone.  Additional  results  have 
also  been  obtained  using  the  modified  Brode  airblast  and  important  comparison  results  will 
be  presented  later.  The  alrblast  loading  which  Impinges  on  the  front  face  is  assumed  to 
act  in  the  gap  cavity  as  well  and  its  action  is  one  that  tends  to  compress  the  closure, 
but  push  the  frame  outward. 

Deformation  Mode — Figure  3-23  shows  the  deformation  mode  of  the  longitudinal  cross- 
section  of  the  structure  with  gap  when  subjected  to  the  AFWL  airblast  and  should  be  com- 
pared with  Figure  3-17,  the  deformed  shape  of  the  structure  without  gap  under  the  same  load- 
ing conditions.  The  snapshots  are  taken  at  10  msec,  after  alrblast  arrival  with  front-on 
incidence  and  for  clarity  the  displacements  are  exaggerated  by  the  same  magnification  factor 
for  both  cases. 

Definite  and  significant  opening  of  the  gap,  as  compared  to  deformation  of  the 
other  portions  of  the  structure,  can  be  readily  observed.  It  is  noted  that  the  opening  is 
larger  at  the  bottom  edge  of  the  closure  than  at  the  upper  edge. 

Another  view  of  the  response  of  the  gap  under  loading  is  given  in  Figure  3-24 
which  shows  the  deformation  portion  of  the  front  face  at  10  msec.  Again,  it  is  clear  that 
the  lower  lip  opens  up  much  more  than  the  upper  lip  and  that  significant  shear  stress  is 
Induced  in  the  headworks  especially  around  the  corners. 

Stresses — to  illustrate  this  point  further,  the  stress  time  history  at  a point  in 
the  lower  corner  of  the  headworks  is  monitored  and  the  result  shown  in  Figure  3-25  in  the 
form  of  generalized  shear  stress  (/j " where  is  the  second  invariant  of  the  deviatoric 
stress  tensor)  versus  mean  normal  stress  (Jj^/3,  where  Jj^  is  the  first  invariant  of  the 
stress  tensor)  plots.  Comparing  the  responses  with  gap  and  without  gap  and  referring  to 
the  approximate  failure  criterion  for  plain  concrete  included  in  the  figure,  it  can  be 
concluded  that  the  presence  of  the  gap  increases  the  possibility  of  failure,  at  least 
locally,  in  the  frame. 

Similar  results  have  been  obtained  for  other  locations  on  the  headworks  front 
face.  It  appears  that  the  four  corner  locations  are  susceptible  to  shear  failure  due  to 
the  presence  of  air  pressure  in  the  gap,  with  the  lower  corners  the  more  vulnerable.  This 
finding  is  consistent  with  the  S-1  post-test  observation  that  "all  four  sides  of  the  head- 
works  (frame)  exhibited  gross  shear  deformation  and  the  top  and  bottom  edges  sheared 
completely  off"  (Reference  11), 

Introducing  the  gap  in  the  model  also  uncouples  the  sides  of  the  closure  from 
the  frame  and  the  total  load  applied  to  the  closure  face  is  now  supported  by  the  bearing 
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ring.  As  shown  In  Figure  3-26  the  peak  bearing  pressures  on  the  upper  and  side  bearing 
ring  are  increased  as  expected.  However,  the  bearing  pressure  on  the  lower  bearing  ring 
shows  a decrease  which  is  rather  surprising.  In  any  event  the  average  peak  bearing  pres- 
sure to  peak  applied  load  ratio  is  below  the  area  ratio  of  2.5.  It  should  be  emphasized 
that  the  pressure  quoted  here  is  actually  the  axial  stress  which  exists  in  the  elements 
which  make  up  the  bearing  ring  and  hence  is  only  an  approximate  measure  of  the  contact 
pressure  in  the  closure/bearing  interface. 

Soil-Structure  Interaction — the  fact  that  the  lower  lip  opens  up  more  than  the 
upper  lip  can  be  due  to  one  or  both  of  the  following  two  reasons.  The  applied  pressure 
is  slightly  higher  in  the  lower  lip  than  the  upper  lip  and  hence  results  in  a slightly 
greater  opening  force.  Also,  the  passive  resistance  to  opening  of  the  lower  lip  due  to 
soil  is  less  than  the  resistance  applied  to  the  upper  lip  due  to  direct  airblast  load. 

The  fact  that  the  same  phenomenon  is  observed  in  the  modified  Brode  case  where  the  upper 
and  lower  lip  (and  front  face)  are  loaded  by  the  same  overpressure  load  suggests  that  the 
latter  is  more  important. 

3.4.  CALCULATION  ID— THREE-DIMENSIONAL  ElJ>iSTIC  CALCULATION  WITH  30°  OBLIQUE  FRONT-ON 
INCIDENCE 

This  calculation  is  performed  to  investigate  the  effect  of  the  angle  of  incidence 
(sometimes  referred  to  as  quartering  attack)  on  structural  response.  The  airblast  approaches 
the  front  face  at  an  angle  of  30°  to  the  longitudinal  plane,  wherein  it  impinges  the  edge  of 
the  front  face  on  the  blast  side  before  sweeping  across  the  face  and  the  exposed  headworks 
surface  and  acting  on  the  back  side. 

Model  and  Airblast  Load — this  is  a full  three-dimensional  calculation  in  which 
no  plane  of  symmetry  exists  and  the  complete  shelter/berm  configuration  needs  to  be  modeled. 
The  head-on  elastic  model  used  in  Calculation  IB  is  expanded  to  include  both  sides  of  the 
soil-structure  about  the  geometric  plane  of  symmetry  and  the  resulting  model  is  illustrated 
in  Figure  3-27. 

For  lack  of  full  three-dimensional  airblast  data,  an  estimate  of  the  loadings 
based  on  HAVE  HOST  test  data  (Reference  11)  and  AFWL  HULL  calculation  data  (Reference  4) 
is  used.  The  longitudinal  component  corresponds  to  that  used  in  the  head-on  calculations 
so  that  any  change  in  shelter  response  observed  in  this  calculation  can  be  attributed  to 
the  change  in  angle-of-lncldence  from  zero  to  30°. 

In  addition  to  the  front  face,  the  blast  side  of  the  headworks  is  subjected  to 
high  overpressure  loading  due  to  reflection  of  the  wave  front  by  the  headworks  and  berm  and 
the  magnitude  assumed  for  this  loading  is  obtained  from  SI  test  data  with  an  approximate 


magnification  factor  of  four.  On  the  back  side  as  well  as  on  top  of  the  headworks  the 
overpressure  drops  drastically  and  the  nominal  free-fleld  pressure  is  assumed  to  act  on 
these  surfaces. 

Deformation  Mode — the  deformation  pattern  of  the  structure  20  msec,  after  first 
arrival  of  the  wave  front  is  given  in  Figure  3-28  for  a cross-section  in  the  vertical 
plane  of  symmetry  (pitch  plane)  and  in  Figure  3-29  for  a cross  section  in  the  horizontal 
plane  of  symmetry  (yaw  plane).  The  pitch  plane  deformation  Is  similar  to  that  observed 
in  the  head-on  calculation.  The  yaw  deformation  consists  of  compression  in  the  longitudinal 
direction  and  lateral  displacement  of  the  headworks  in  the  direction  of  wave  propagation. 

Again,  the  displacements  are  magnified  to  illustrate  the  deformation  mode. 

( 

( Stresses — axial  stress-time  histories  (Figures  3-30  and  3-31  for  the  pitch  plane) 

are  similar  to  those  obtained  in  the  head-on  calculation.  The  effect  of  quartering  shows 
up  mainly  in  stresses  and  bending  moments  in  the  yaw  plane.  As  shown  in  Figures  3-32  and 
3-33,  in  addition  to  exhibiting  the  expected  amplification  as  it  progresses  along  the  struc- 
ture, the  axial  stress-time  history  in  the  yaw  plane  also  shows  the  effect  of  arrival  time 
since  the  wave  front  impinges  first  on  the  blast  edge  then  progressively  sweeps  across  the 
closure  and  front  face  before  the  back  edge  experiences  the  loading.  This  imbalance  in 
axial  stress  in  the  yaw  plane  due  to  wave  front  arrival  times  generates  a yaw  bending 
moment  in  such  a way  that  the  headworks  is  initially  torqued  toward  the  blast.  Later  on 
the  lateral  load  from  over-pressure  acting  on  the  blast  side  of  the  headworks  tend  to 
counter  balance  this  torque  and  turns  the  headworks  with  the  wind.  This  behavior  is  shown 
in  Figure  3-34,  the  bending  moment  being  measured  by  the  difference  between  the  axial 
stresses  on  the  blast  and  back  sides.  The  maximum  yaw  bending  moment  obtained  is  about 
1.7  X 10*  in. -lbs.,  approximately  half  of  that  observed  in  the  pitch  plane. 

The  effect  of  quartering  on  shear  stress  distribution  as  measured  by  /j  in  the 
headworks  and  shelter  is  also  investigated  and  found  not  to  be  significant  in  general,  A 
typical  comparison  of  the  head-on  calculation  and  quartering  calculation  with  respect  to 
shear  is  shown  in  Figure  3-35.  Each  data  point  in  the  figure  corresponds  to  the  maximum 
shear  stress  measured  relative  to  the  reference  failure  criterion  at  a location  along  the 
structure.  In  both  the  head-on  and  the  quartering  cases  shear  stress  is  amplified  as  one 
moves  along  the  longitudinal  direction  of  the  shelter  in  much  the  same  way  that  axial 
stress  is  amplified  due  to  the  decrease  in  load  bearing  area.  The  quartering  calculation 
gives  slightly  higher  shear  stresses  than  those  obtained  from  the  head-on  calculation;  the 
magnitude  of  the  difference  is  however  not  significant. 
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NONDIMENSIONALIZED  OVERPRESSURE 


Figure  3-5a.  Deformation  mode  T = 40  msec.,  side-on  incidence,  Calculation  lA 


Figure  3-5b.  Deformation  mode  at  T = 60  msec.,  side-on  incidence,  Calculation  lA 


Fourier  spectrum  of  de forma tional  mode  at  T = 40  msec.,  slde-on  incidence,  Calculation  lA 
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Figure 
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Figure  3-8. 


Figure  3-9. 


ontal  and  vertical  stress/time  histories  in  soil,  rear  springline,  side-on  incidence 
latlon  lA 


Figure  3-11.  Three-dimensional  finite 
Calculations  IB  and  1C 


Figure  3-12a.  Cross-sectional  view  of  three-dimensional  model,  front  face  and  headworks 


Figure  3-14.  Comparison  of  peak  overpressure  profiles,  frcnt-on  incidence.  Calculations  IB  and  1C 


Figure  3-15b.  Comparison  of  overpressure  waveforms,  top,  front-on  inci 


Front  Face  Exposed 
to  Airblast 


1-16.  Comparison  of  peak  overpressure  distribution  on  front  face 


Figure  3-18.  Axial  st ress-t ime  histories  along  shelter  length,  upper  headworks,  front-on  incidence 
Calculation  IB 


APPROX.  FAILURE  CRITERIA 
FOR  PLAIN  CONCRETE 


Plot  of  first  and  second  stress  Invariants  at  front  face,  gap,  front-on  incidence 
Calculation  1C 


APPROX.  FAILURE  CRITERIA 
FOR  PLAIN  CONCRETE 


front-on  incidence.  Calculation  1C 


Cal 


Figure  3-28.  Deformation  mode,  vertical  plane  of  symmetry,  30  oblique  front-on  Incidence 
Calculation  ID 
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Figure  3-30.  Variation  in  axial  stress  along  top  of  structure,  30°  oblique 
front-on  incidence,  Calculation  ID 
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Figure  3-31.  Variation  in  axial  stress  along  bottom  of  structure,  30°  oblique 
front-on  incidence,  Calculation  ID 
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Figure  3-33.  Variation  in  axial  stress  along  springline,  back  side  of  structure, 
30°  oblique  front-on  incidence.  Calculation  lA 


Incidence,  Calculation  ID 


Flgur;2  3-35.  Comparison  of  shear  stress  distribution  In  longitudinal  direction,  upper  shelter,  30°  oblique 
front-on  incidence.  Calculation  ID 


TABLE  3-1 


SHELTER  "RING"  IN  VACUO  FREQUENCIES  (HZ) 


Flexural  Mode 

Finite  Element  Model 

Thin  Shell  Theory 

n = 2 

25-27. 

26 

n = 3 

72. 

73. 

II 

c 

131-133. 

141. 

n = 5 

202. 

227. 

n = 6 

280. 

333. 

n = 7 

359. 

459. 

Breathing  Mode 

o 

11 

c 

138. 

134. 
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TABLE  3-2 


ELASTIC  PROPERTIES  OF  SOIL  AND  STRUCTURE 

C 

E P 


Material 

(psf) 

V 

(fps) 

Native  Soil 

1.016E7 

.3 

2,000 

Berm 

2.46E6 

.1 

900 

Backfill 

2.46E6 

.1 

900 

Concrete 

6.48E8 

.2 

10,000 

SECTION  IV 

INTERPRETATION  OF  PHASE  1 RESULTS 

Major  results  of  Calculations  lA  through  ID  all  based  on  the  SI  configuration  are 
discussed  herein.  The  discussion  begins  with  a brief  summary  of  the  structural  response 
mode  in  order  to  lead  into  the  conclusions  which  follow.  The  response  of  the  SI  configura- 
tion is  used  only  to  Illustrate  certain  observations  on  the  methods  of  the  dynamic  analysis 
of  a shelter-like  structure.  Note  also  that  the  Phase  1 results  are  obtained  through  elas- 
tic calculations  and  the  conclusions  are  drat.m  based  on  the  elastic  assumption.  Some  con- 
clusions are  modified  as  inelastic  calculation  results  from  Phase  2 become  available. 

4.1.  DEFORMATION  MODE 

The  main  deformation  mode  of  the  shelter  under  head-on  loading  is  similar  to  that 
of  a beam  column  in  an  elastic  medium  (see  for  example  Figure  3-17).  The  tremendous  load 
acting  on  the  front  face  dominates  the  response  of  the  structure  at  early  times  and  produces 
two  major  effects:  high  axial  stresses  in  the  headworks  and  particularly  the  tube  (Figures 
3-18  and  3-19)  and  bending  of  the  headworks  as  a result  of  the  offset  between  the  center  of 
pressure  and  the  center  of  stiffness  (Figure  3-20).  The  alrblast  on  the  top  of  the  berm 
tends  to  push  the  headworks  downward  and  contributes  to  the  bending  effect.  The  effect  of 
the  overburden  load  lags  the  effect  of  the  front  load  because  of  the  necessary  transmission 
time  through  the  backfill. 

When  the  airblast  incidence  is  oblique,  there  exists  in  addition  to  the  deforma- 
tion described  above  a yaw  deformation  consisting  of  compression  in  the  longitudinal  direc- 
tion and  lateral  displacement  (Figure  3-29)  and  bending  of  the  headworks  in  the  direction 
of  wave  propagation.  The  effect  of  the  oblique  incidence  is  first  manifested  as  the  effect 
of  alrblast  arrival  time  at  the  front  face,  torquing  the  headworks  initially  into  the  blast. 
Later  on  the  lateral  load  acting  on  the  blast  side  of  the  headworks  tends  to  counterbalance 
this  torque  and  turns  the  headworks  with  the  blast  (Figure  3-34). 

Tlie  longitudinal  and  lateral  bendings  of  the  headworks  are  transmitted  to  the 
tubular  section.  Consequently  the  tube  section  closest  to  the  headworks  is  under  severe 
distress  especially  at  the  crown.  In  addition,  although  only  substantiated  by  a two- 
dimensional  calculation  the  tube  undergoes  In-plane  deformations  due  to  the  overburden 
load.  This  in-plane  deformation  Is  dominated  by  the  breathing  and  oval Ing  modes  (Figure  3-5). 

Other  possible  areas  of  distress  have  been  identified  as  the  Inside  corners  of 
the  frame;  the  distress  is  local  and  is  caused  exclusively  by  the  air  pressure  in  the  gap, 
I.e.,  the  space  between  the  closure  and  the  frame  (see  Figure  3-25). 
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The  bending  moment  at  the  center  of  the  closure  Is  compared  to  the  values  given 
by  plate  theory  In  Figure  4-1.  It  can  be  summarized  that  the  closure  In  Its  elastic  state 
behaves  like  a thick  plate  In  flexure  with  fixed  edges  when  the  gap  Is  not  present  and  Is 
simply-supported  otherwise. 

The  velocity-time  histories  for  a point  s-t  the  center  of  the  closure  in  the  mono- 
lithic front  face  (i.e..  no  gap)  are  shown  in  Figure  4-2  and  contain  a strong  component 
with  a period  of  about  8 msec.  The  maximum  midspan  deflection  relative  to  the  support  is 
less  than  an  inch.  Corresponding  time  histories  for  the  case  with  gap  are  not  available 
since  that  particular  calculation  is  terminated  at  10  msec,  after  alrblast  first  arrival. 

4.2.  COMPUTER  MODELING  CONSIDERATIONS 

It  follows  from  the  foregoing  discussions  that  several  regions  of  the  shelter 
need  to  be  modeled  in  considerable  detail.  These  are:  the  transition  region  between  the 
headworks  and  the  tube  including  the  front  portion  of  the  tube,  the  gap  between  the  closure 
and  frame,  the  inside  corners  in  the  frame  and  the  front  face  and  closure  itself. 

Details  in  the  transition  region  and  tube  should  be  sufficient  to  model  the 
longitudinal  stress  transfer/amplification  and  the  severe  distress  induced  by  bending  of 
the  headworks.  The  gap  between  the  closure  and  frame  should  be  included  if  air  pressure 
is  allowed  to  get  into  this  cavity,  i.e.,  no  gap  seal  is  anticipated.  The  mesh  in  upper 
and  lower  corners  on  the  inner  surface  of  the  frame  should  be  sufficiently  fine  to  reflect 
the  high  stress  gradient  which  exists.  Finally  the  front  face  details  should  adequately 
represent  the  spatial  variations  of  the  overpressure  distribution. 

It  has  also  been  found  during  the  course  of  the  study  that  the  soil  region 
immediately  in  front  of  the  shelter  plays  an  Important  role  in  the  shelter  response.  The 
alrblast  load  acting  on  this  surface  changes  very  rapidly  with  distance  from  the  front 
face  as  the  reflected  overpressure  decays  rapidly  from  its  full  value  at  the  front  face. 
This  distribution  can  only  be  reproduced  with  a fine  mesh  for  the  soil  surface  in  question. 

Of  course  a fine  mesh  is  also  necessary  to  allow  accurate  transmission  of  the  load  into 

the  buried  portion  of  the  headworks. 

A fine  surface  subdivision  may  be  necessary  for  another  reason.  It  is  common 

practice  in  most  codes  to  distribute  the  load  acting  on  a surface  equally  to  nodes  which 

define  that  surface.  Consider  the  discretized  soil  surface  adjacent  to  the  front  face 
which  has  nodes  on  the  structure  and  nodes  on  the  soil;  it  this  surface  is  coarse  (large) 
part  of  the  high  soil  surface  load  will  act  (erroneously)  on  the  structure,  trying  to 
bring  it  downward  together  with  the  soil.  This  pitfall  is  to  be  avoided,  particularly  in 
an  elastic  calculation  where  the  strength  of  the  soil  material  is  not  limited. 
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Another  Important  aspect  of  the  model  is  the  modeling  of  the  material  properties. 
For  example,  in  an  elastic  analysis,  material  modeling  Involves  choosing  the  elastic  material 
constants  (bulk  and  shear  moduli)  which  best  approximate  the  real  material  behavior  in 
particular  situations.  Although  it  is  not  considered  as  part  of  the  scope  of  the  present 
study,  the  success  of  the  methodology  established  here  for  the  dynamic  analysis  of  shelter- 
like structures  depends  strongly  on  adequate  material  modeling  for  the  reinforced  concrete, 
steel  and  the  soil  medium  as  is  true  for  any  analysis  of  this  kind.  In  particular  as  will 
be  shown  in  Section  VI,  the  effectiveness  of  the  methodology  In  predicting  closure  response 
depends  on  an  accurate  dynamic  model  representation  of  the  closure  materials. 

4.3.  AIRBLAST  CHARACTERISTICS  AND  SHELTER  RESPONSE 

As  far  as  their  effects  on  shelter  response  is  concerned,  the  characteristics  of 
the  overpressure  wave  form  can  be  discussed  under  two  separate  groups:  spatial  and  temporal. 
Spatial  characteristics  refer  to  the  distribution  of  peak  overpressures  acting  on  the  struc- 
ture and  ground  surface.  Temporal  characteristics  refer  to  the  details  of  the  time  history 
and  include  parameters  such  as  arrival  time,  rise  time,  dominant  frequency  and  overpressure 
impulse.  These  two  sets  of  characteristics  are  interrelated  and  it  is  often  difficult  to 
isolate  the  effects  of  one  from  the  other.  The  temporal  parameter  found  to  be  significant 
is  the  arrival  time  of  the  peak  overpressure  at  a point  and  equivalently  the  relative  timing 
of  the  various  peak  load  arrival  times.  Spatial  characteristics  are  found  to  be  important 
in  some  situations  and  not  so  in  others. 

Before  going  into  a discussion  of  these  characteristics,  however,  it  is  emphasized 
that  the  results  obtained  in  this  respect  are  limited  because  of  the  lack  of  authentic  air- 
blast  data.  What  little  comparison  has  been  done  is  between  the  AFWL  19.8041  (Reference  4), 
19.8042  (Reference  5)  loads  obtained  under  two-dimensional  conditions  and  the  modified 
Brode.  Accordingly,  this  aspect  of  the  study  is  limited  in  scope  and  the  conclusions  tenta- 
tive. Furthermore,  the  conclusions  are  based  on  elastic  results  which  tend  to  exaggerate 
the  response  sensitivity. 

The  elastic  response  of  the  shelter  depends  a great  deal  on  the  peak  applied 
load  and  its  distribution  on  the  loaded  surface  area.  This  is  so  because  of  the  peculiar 
geometry  of  the  shelter,  with  its  exposed  front  face  and  sides,  top,  gap,  berm,  overburden, 
etc. , and  the  sensitivity  of  the  peak  overpressure  to  this  geometry,  the  amplification  fac- 
tor being  such  a nonlinear  and  sensitive  function  of  the  reflection  and  diffraction  angles. 
This  sensitivity  is  best  illustrated  by  comparing  the  peak  overpressure  distributions  used 
in  the  two-dimensional  studies  (Figure  3-3)  and  three-dimensional  studies  (Figure  3-14). 

In  the  three-dimensional  studies,  the  peak  overpressure  suffers  a big  jump  in  value  as  the 


airblast  encounters  the  door;  the  factor  of  reflection  ranges  from  7 to  9.  This  peak  then 
decreases  rapidly  as  the  airhlast  rounds  the  top  and  side  edges  of  the  structure  and  finally 
regains  the  free  air  profile  some  distance  down  the  herm.  Meanwhile,  the  reflected  airhlast 
acts  on  the  soil  medium  in  front  of  the  shelter  as  well  with  peak  value  decaying  rapidly 
with  distance  from  the  front  face. 

To  the  structure  the  front  load  collected  at  the  front  face  must  he  supported  hy 
the  headworks  and  the  medium  surrouncing  it.  High  axial  stresses  result  and  the  situation 
is  made  more  acute  as  the  stress  wave  propagates  along  the  structure  due  to  a drop  in  load 
hearing  area  in  going  from  the  headworks  to  the  tube.  This  effect  is,  without  doubt,  the 
single-most  dominant  effect  of  airblast  characteristics.  Hence,  the  front-face  reflection 
coefficient  of  the  surface  overpressure  is  an  important  airblast  parameter. 

The  front-face  reflection  coefficient  is  not  the  only  important  parameter  related 
to  the  front  face.  Its  distribution  about  the  front  face  is  another  significant  factor. 

To  Illustrate  this  point,  the  modified  Erode  distribution  (see  Figure  3-16)  which  has  the 
same  high  peak  overpressure  (5,400  psi)  as  the  AFIVL  8041  load,  but  is  uniformly  applied 
over  the  entire  front  face  is  used  to  generate  results  for  comparison.  This  distribution 
is  different  enough  from  the  AFIJL  19.8041  distribution  to  yield  some  interesting  response 
sensitivity  observations.  Figure  3-16  shows  a portion  of  the  front  face  shielded  by  soil 
and  hence  not  subjected  to  direct  airblast  loading.  On  that  portion  which  is  exposed  is 
plotted  the  peak  overpressure  distributions.  The  dotted  line  represents  the  AFVL  19.8041 
pressure  and  the  dashed  line  represents  the  modified  Erode  pressure  constant  over  the  front 
face.  The  volume  under  each  surface  is  the  total  axial  load  acting  on  the  front  face  in  each 
case  and  the  moment  about  the  center  of  the  front  face  can  be  measured  accordingly. 

Due  to  the  presence  of  the  layer  of  soil  in  the  lower  front-face  area  the  center 
of  pressure  is  offset  from  the  shelter  center  line.  This  offset  in  the  modified  Erode 
description  is  larger  than  that  in  the  AFWL  description  because  of  the  uniform  assumption. 

So  is  the  resultant  load.  This  offset  of  about  4 feet  in  combination  with  the  higher  air- 
blast  load  (4.7  10°  lbs.  for  uniform  distribution)  generates  a moment  of  2.4  x lo’°  in.- 

Ibs. , which  turns  the  headworks  clockwise.  Hence  it  is  not  surprising  that  the  deformation 
mode  of  the  structure  (Illustrated  in  Figure  4-3)  due  to  the  modified  Erode  load  shows  both 
a greater  axial  compression  and  higher  bending  of  the  headworks.  In  fact,  the  bending  ac- 
tion is  so  severe  that  there  is  a lift-off  of  the  headworks  at  the  toe,  a phenomenon  also 
observed  In  the  SI  test  (Reference  11).  The  deformation  mode  of  the  same  structure  and  at 
the  same  time  for  the  AFWL  load  is  also  included  for  comparison.  It  also  follows  that 
compared  to  the  AFWL  load  the  modified  Erode  load  should  also  produce  higher  mean  axial 
compression  (14  ksl  versus  7 ksi)  and  higher  bending  moment  in  the  tube  (1.3  x 10* ° in. -lbs. 
versus  0.36  x 10*°  in. -lbs.)  as  shown  in  Figure  4-4.  This  load/response  sensitivity 

4-4 


Is  the  most  evident  because  the  load  is  not  only  directly  applied  to  the  structure,  but 
also  the  most  dominant. 
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A similar  conclusion  can  be  drawn  regarding  the  load  distribution  on  the  exposed 
side  of  the  headworks  in  the  slde-on  incidence  case,  although  the  actual  calculation  has 
not  been  made.  An  indication  that  this  should  follow  the  trend  of  the  axial  load  in  the 
head-on  incidence  is  given  by  results  of  the  oblique  front-on  calculation.  The  relevant 
parameters  are  now  the  lateral  bending  and  shear  in  planes  perpendicular  to  the  longitudinal 
axis. 

When  the  loading  is  not  direct,  however,  this  sensitivity  is  greatly  diminished. 

Such  are  the  effects  of  the  overburden  load  on  the  headworks  and  tube.  First  of  all  the 
overburden  load  effect  is  delayed  somewhat  because  of  the  berm.  For  instance,  for  the 
material  assumed,  it  takes  a little  more  than  10  msec,  for  the  airblast-lnduced  ground 
motion  to  traverse  the  overburden  and  reach  the  tubular  section  of  the  shelter.  Accordingly, 
the  magnitude  of  this  loading,  small  to  begin  with  when  compared  to  the  front  load,  is  further 
attenuated  due  to  propagation  through  the  backfill  material.  An  Indication  of  this  decrease 
in  sensitivity  can  be  obtained  by  looking  at  the  results  of  the  elastic  two-dimensional 
calculation. 

As  shown  in  Figure  3-3,  the  peak  overpressure  at  the  blast  side  of  the  berm  ex- 
periences amplification  of  a factor  of  two  above  the  pressure  for  flat  ground  (Erode).  This 
load,  however,  must  travel  through  about  80  feet  of  backfill  material  in  order  to  reach  the 
structure  and  hence,  as  one  would  expect,  the  effect  of  this  amplification  of  two  on  struc- 
tural response  is  small.  For  instance,  the  tube  response  shows  slightly  higher  hoop  com- 
pressions in  the  structure,  but  not  in  the  same  proportion  as  the  amplification  (see 
Figure  4-5).  The  bending  stress  due  to  the  ovaling  mode  is  also  quite  insensitive  to  the 
details  of  the  overpressure  distribution  also  (3,800  psi  and  2.5  x lO®  in. -lbs. /in.  for 
the  AFWL  load  and  3,200  psi  and  2.6  x 10®  in. -lbs. /in.  for  the  Erode  load).  In  particular, 
when  more  realistic  soil  properties  are  used  to  include  hysteresis  and  damping,  any  sensi- 
tivity shown  in  the  elastic  calculations  will  diminish  even  further. 

With  regard  to  temporal  characteristics,  there  are  at  least  two  important  para- 
meters which  need  to  be  addressed:  the  rise/decay  time  associated  with  the  overpressure 
peak  and  the  Impulse  associated  with  the  load.  From  a frequency  viewpoint,  the  former  is 
considered  a high  frequency  effect  and  the  latter  belongs  to  the  medium/low  frequency  regime. 

The  importance  of  these  parameters  and  the  sensitivity  of  the  structural  response 
to  them  again  depends  on  the  loading  configuration.  In  the  front-on  case,  the  direct  applica- 
tion of  the  load  to  the  exposed  surface  of  the  structure  affords  the  best  conditions  for 
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the  transfer  of  the  high  frequency  vibration  to  the  structure.  Potential  problems  are  j 

resonance  of  the  closure,  if  it  remains  elastic,  and  shock  isolation  of  mlsslle/equipment 
inside  the  headworks.  Farther  down  the  shelter,  the  pulse  wave  form  is  dispersed  and  dis- 
torted due  to  the  geometry  of  the  structure  and  the  medium;  the  frequency  content  is  modi-  i 

fled.  When  loading  is  not  direct,  the  backfill  and  soil  mass  surrounding  the  structure 
will  filter  out  high  frequency  contents  of  the  overpressure  and  structural  response  is 
expected  to  be  much  less  sensitive  than  in  the  direct  loading  case. 

I [ 

The  overpressure  Impulse  appears  to  be  the  most  difficult  parameter  to  be  Isolated 
and  studied.  It  is  usually  a long  time  effect  with  Influence  on  the  maximum  displacement, 
permanent  set  and  ductility  ratio,  etc.  An  analysis  of  this  parameter  is  made  more  difficult 
because  the  results  obtained  from  such  a calculation  are  usually  masked  by  different  loads 
originating  from  different  stations.  It  is  probably  more  pertinent  to  address  this  para- 
meter in  an  inelastic  calculation,  as  is  done  in  Phase  2 in  discussing  the  Inelastic 
response  of  the  closure. 
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Figure  A-1.  Bending  moment  coefficients  at  center  of  closure,  front-on  Incidence,  comparison  of  finite 
element  and  plate  theory  results 
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Figure  4-4a.  Axial  stress  time  histories,  tubular  portion  behind  headworks,  modified  Erode,  front-on 
incidence 


Figure  A-Ab.  Axial  stress  time  histories,  tubular  portion  half  tube  diameter  behind  headworks,  modified 
Erode,  front-on  incidence 


SECTION  V 


PHASE  2 RESULTS 

Since  the  Initiation  of  this  study,  the  design  of  the  MAP  shelter  prototype  has 
undergone  several  major  changes;  it  was  decided  that  the  Phase  2 studies  should  parallel 
the  concurrent  design  of  the  shelter,  designated  the  S4.  Furthermore,  the  scope  of  the 
Phase  2 studies  was  modified  In  order  to  provide  analytic  support  to  the  testing  of  a 
one-half  size  shelter  in  the  HAVE  HOST  DABS  (Reference  12)  series.  This  test  was  even- 
tually cancelled  and  replaced  by  the  HEST  S4  test  (Reference  13).  The  structure,  however. 
Is  basically  the  same  one  and  will  be  referred  to  simply  as  the  S4  structure  hereinafter. 
The  full  size  S4  structure  is  studied  in  Phase  2. 

Some  of  the  major  changes  in  the  design  from  SI  to  S4  are  (see  Figures  5-1  and 

5-2): 

• A smaller  front  face  which  slants  at  10°  to  the  vertical. 

• A smaller  shelter  diameter  and  a thinner  wall. 

• A different  closure/frame  design;  the  closure  is  a thick  composite  slab  (con- 
crete supported  by  steel  side  and  back  plates)  with  bevelled  edges. 

• The  bearing  ring  (closure  support)  area  is  smaller. 

• The  headworks  remains  basically  rectangular,  but  has  a thicker  wall  on  the 
hinge  side  to  accommodate  the  power  actuator. 

• The  center  line  of  the  rectangular  headwork  cavity  is  offset  from  the  center 
line  of  the  tubular  portion  of  the  shelter. 

• Only  the  front  face  and  a very  small  portion  of  the  top  of  the  headworks  are 
exposed.  The  top  is  covered  by  overburden  and  the  sides  by  two  wing  walls  and  fill. 

These  changes  and  the  requirement  to  include  as  much  detail  as  possible  of  the 
closure/headworks  resulted  in  a formidable  finite  element  model.  The  advantage  of  such  a 
model  Is  that  detailed  response  of  the  structure,  especially  the  closure/headworks  sec- 
tion can  be  obtained;  the  disadvantage  is  that  small  integration  time  steps  must  be  used. 
Even  with  the  subcycling  capability  of  TRANAL,  this  requirement  still  results  in  long 
computation/ run  time. 

All  Phase  2 calculations  use  basically  the  same  model  which  Includes  the  gap 
between  closure  and  frame.  They  are  all  made  in  the  head-on  loading  condition  with  the 
same  overpressure  description.  Major  results  of  these  calculations  will  be  described  in 
this  section;  a detailed  interpretation  of  the  Phase  2 findings  is  given  in  Section  VI. 
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5.1.  CALCULATION  2A— THREE-DIMENSIONAL  INELASTIC  CALCULATION  WITH  HEAD-ON  INCIDENCE 
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Model  and  Alrblast  Load — symmetry  Is  assumed  about  the  center  vertical  plane  of 
the  tube  and  the  model  includes  only  half  (hinge  side)  of  the^shelter/berm  configuration. 
The  finite  element  model  which  consists  of  about  18,000  hexahedrons  is  shown  in  Figure  5-3. 
As  can  be  seen  from  this  figure  special  attention  has  been  given  to  details  in  the  closure, 
headworks  and  transition  region.  In  addition  to  the  gap,  the  steel  side  and  back  plates 
encasing  the  concrete  closure  are  modeled  using  two  elements  across  the  plate  thickness. 

The  steel  bearing  ring  and  the  steel  plates  in  the  frame  are  modeled  likewise.  The  clo- 
sure Is  assumed  rigidly  bonded  to  the  bearing  ring  where  contact  is  made.  Two  elements 
are  used  across  the  thickness  of  the  shelter  cylindrical  wall. 

The  In  situ  soil  and  backfill  are  modeled  as  "CAP"  materials;  the  corresponding 
CAP  parameters  are  given  in  Table  5-1.  The  in  situ  soil  model  as  presented  is  identical 
to  that  designated  as  DRY  SAND-1  in  Weldlinger  Associates'  ground  motion  study  DNA  001-77- 
C-0036  (Reference  9).  DRY  SAND-1  has  a density  of  110  pcf  and  an  uniaxial  loading  modulus 
corresponding  to  a wave  speed  of  1,500  fps  and  an  unloading  wave  speed  of  3,600  fps.  Its 
uniaxial  behavior  is  illustrated  in  Figure  5-A. 

The  backfill  model  is  based  on  data  obtained  from  the  AFWL  (Reference  7)  and  cor- 
responds to  the  so-called  medium  backfill  of  the  AFWL  data.  The  uniaxial  behavior  of  the 
model  is  compared  in  Figure  5-4b  with  the  AFWL  data;  the  loading  wave  speed  is  800  fps  and 
the  unloading  wave  speed  3,600  fps. 

The  concrete  and  steel  are  modeled  as  elastic  perfectly  plastic  materials,  the 
concrete  having  an  exponential  yield  surface  based  on  a compressive  strength  of  6,000  psl 
and  the  steel  a von  MIses'  yield  surface  corresponding  to  A36  steel.  The  model  parameters 
are  listed  in  Table  5-1.  In  the  concrete  model  the  tension  cut-off  point  is  varied  based 
approximately  on  the  amount  of  reinforcing  steel  that  exists  in  different  parts  of  the 
prototype  structure  (i.e.,  headworks,  frame,  closure  and  tube)  and  the  basic  yield  surface 
is  kept  unchanged. 

Since  detailed  airblast  loading  on  the  S4  shelter  configuration  does  not  exist, 
a combination  of  test  and  handbook  data  is  used  for  this  calculation  based  on  recommenda- 
tions from  AFWL  (Reference  8).  A sample  of  these  data  are  given  in  Figure  5-5  where  they 
are  compared  with  the  axlsvmmet rlc  AFWL  HULL  results.  The  loading  at  the  front  face  is 
found  Identical  to  the  HULL  values  and  close  to  that  measured  in  the  HAVE  HOST  SI  test. 

The  integration  time  step  used  in  TRANAL  can  be  different  for  different  groups  of 
elements  called  zones.  Within  each  zone,  however,  the  time  step  Is  constrained  by  the 
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minimum  transit  time  across  an  element  in  that  zone,  In  Calculation  2A  the  time  step  for 
the  zone  containing  the  closure  details  is  necessarily  small  due  to  the  small  size  (2.5 
inches)  of  the  elements  used  to  model  the  steel  side  and  back  plates.  Consequently  the 
computation  time  per  major  cycle  is  long:  about  50  CDC  7600  CPU  seconds  for  0.1  msec, 
of  real  time  and  approximately  18,000  elements.  (This  figure,  of  course,  would  be  higher 
if  not  for  the  subcycling  capability  of  the  code.)  Because  of  the  formidable  run  time 
only  about  10  msec,  of  structural  response  time  has  been  obtained  for  Calculation  2A. 
Although  not  long  enough  to  fully  examine  the  behavior  of  the  tube  portion  of  the  struc- 
ture this  response  time  is  sufficient  to  yield  results  for  an  analysis  of  the  headworks 
and  closure. 

Deformation  Mode — an  illustration  of  the  deform.irion  pattern  of  the  headworks 
is  given  in  Figure  5-6  with  magnified  displacements  for  points  on  the  longitudinal  plane. 

It  is  seen  that  the  early  time  response  of  the  structure  is  basically  longitudinal  compres- 
sion. Front  loading  acting  at  the  closure  tends  also  to  push  it  downward  In  addition  to 
backward  due  to  the  10°  incline  and  the  bearing  support  is  forced  to  move  with  it  since 
perfect  bond  is  assumed  in  the  model.  The  overall  downward  movement  elsewhere  in  the 
headworks,  however,  is  not  significant.  This  is  because  at  early  times  the  soft  over- 
burden effectively  isolates  the  structure  from  top  airblast  loads. 

In  Figure  5-7  are  shown  the  displaced  profiles  of  the  b.ukp 1 at e vt  rlical  center 
line  at  different  times.  Profiles  for  the  horizontal  center  line  are  given  in  Figure  5-8. 
The  deformation  pattern  is  one  of  punching,  or  shear,  where  the  plate  yields  along  the 
support  edges,  but  retains  a relatively  uniform  (flat  plate)  displacement  elsewhere. 

The  development  of  this  deformation  mode  is  illustrated  in  Figure  5-9  which  shows 
the  growth  of  inelastic  regions  with  time  in  the  hack  plate.  It  is  apparent  that  the  area 
near  the  bottom  support  yields  first,  followed  by  areas  near  the  side  support  and  finally 
those  near  the  top  support. 

The  deformation  pattern  In  the  concrete  portion  of  the  composite  closure  is 
shown  in  Figure  5-10  where  three  patterns  arc  given  for  each  action  time  corresponding  to 
(from  left  to  right)  the  top,  mid  and  bottom  layers  which  make  up  the  concrete  slab  model. 
The  bottom  layer  undergoes  inel.istic  deformation  first  which  is  then  seen  to  propagate  to 
the  mid  and  top  layers  at  later  times.  It  is  also  clear  that  inelastic  deformation  is 
initiated  In  the  concrete  before  it  Is  Initiated  in  the  back  plate,  giving  some  evidence 
to  the  theory  that  shear  deformation  propagates  from  hack  to  .front  of  the  concrete  slab; 
it  is  only  after  the  completion  of  the  so-called  shear  crack  that  the  load  Is  transferred 
to  the  backplate  causing  it  to  yield  also.  Results  obtained  appear  to  be  consistent  with 
this  theory  although  the  deformation  pattern  is  more  complex  due  to  dynamic  and  wave 
propagation  effects. 
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Velocity-time  history  for  the  closure  center  of  gravity  (eg)  is  shown  in  Figure 
5-11.  It  is  a predominantly  longitudinal  motion  with  some  downward  component,  due  to  the 
10°  incline  of  the  closure. 

Stresses — the  stress  Invariant  plot  of  Figure  5-12  is  typical  of  result  obtained 
for  points  in  the  concrete  portion  of  the  closure  where  it  becomes  Inelastic.  Notice  that 
the  mean  pressure  (Jj^/3)  remains  compressive  and  the  inelastic  deformation  is  due  to  shear. 
For  regions  in  the  closure  which  remain  elastic  the  stress  Invariant  plot  looks  like  Fig- 
ure 5-13. 

Stress  pattern  in  the  backplate  falls  into  two  categories,  supported  or  unsupport- 
ed. The  supported  or  edge  region  behavior  is  shown  in  Figure  5-14;  the  unsupported  regions 
is  dominated  by  the  membrane  tension  and  the  resulting  stress  invariant  plot  has  the  char- 
acteristics of  Figure  5-15.  The  maximum  shear  allowed  (/T^)  is  20,800  psl  and  corresponds 
to  a uniaxial  yield  stress  of  36,000  psl  for  A36  steel. 

Experience  from  Phase  1 results  indicates  that  other  vulnerable  areas  in  the  struc- 
ture are  the  inner  comers  of  the  frame.  Figures  5-16  and  5-17  are  the  stress  invariant 
curves  for  the  upper  and  lower  comers,  respectively,  of  the  steel  ring  lining  the  frame. 

Two  elements  are  used  to  model  the  ring  thickness,  and  in  both  the  upper  and  lower  corners, 
the  inner  ring  element  (which  is  in  direct  contact  with  the  gap  pressure)  sustains  hoop 
tension  while  the  outer  ring  element  (which  is  adjacent  to  the  frame  concrete)  sustains 
compressive  stresses  as  the  ring  is  bent  and  the  corner  "rounded."  Yielding  is  observed 
in  both  cases. 

Minor  inelastic  deformation  is  also  observed  in  the  concrete  portion  of  the  frame 
near  the  corners  and  center  gussets  as  Figures  5-18  through  5-20  indicate. 

The  stress  Invariant  plot  of  Figure  5-21  is  typical  of  the  results  obtained  in 
the  transition  section  which  remains  elastic,  due  to  the  beneficial  confining  effect  of 
the  neighboring  soil.  As  the  stress  wave  travels  onward  to  the  tubular  section,  the  de- 
crease in  load  bearing  area  results  in  an  Increase  in  the  longitudinal  stress.  This  axial 
compression  is  sufficient  to  cause  inelastic  deformation  in  the  tubular  section  behind 
the  transition  section.  Figure  5-22  is  typical  of  the  results  obtained  at  the  crown  and 
springline.  The  Invert  portion  remains  elastic  at  this  early  time. 

This  amplification  of  the  longitudinal  stress  as  shown  in  Figures  5-23  and  5-24 
has  been  mentioned  in  previous  discussions  (Sections  3.2  and  3.4).  In  the  case  of  the 
Inelastic  calculation,  the  amplification  factor  measured  is  less  than  2 which  is  much 
less  than  the  ratio  of  the  load  bearing  areas  (about  6 for  the  S4  configuration).  There 
are  two  possible  explanations:  the  axial  stress  which  can  be  sustained  in  the  tubular 


section  Is  limited  for  an  inelastic  material  depending  on  its  yield  strength  and  the 
lateral  confinement.  Secondly,  the  front  load  which  acts  on  the  closure  is  not  completely 
transmitted  to  the  transition  section  due  to  the  shear  or  punching  deformation  at  the  sup-  ^ 

port  edges.  More  insight  into  the  former  can  be  gained  by  looking  at  results  of  Calcula-  J 

tlons  2B  and  2C  to  be  presented  later  (Sections  5.1  and  5.2).  The  matter  of  the  bearing  ! 

1 

load  is  examined  in  the  next  paragraph.  | 

The  ratio  of  th^  closure  area  versus  bearing  area  is  approximately  12.  Hence,  I 

it  follows  that  in  static  uniform  loading  of  the  closure,  the  average  bearing  pressure 

is  12  times  that  of  the  loading,  or  a peak  of  64,000  psi  in  this  case.  The  bearing  pres- 
sure obtained  in  Calculation  2A  is  shown  in  Figure  5-25  for  points  located  on  the  top,  ‘ 

bottom  and  side  of  the  bearing  ring.  The  peak  pressure  is  higher  at  the  top  and  bottom 
mainly  because  the  bearing  area  is  less  there  than  on  the  side.  In  any  event,  its  value  1 

of  30,000  psi  is  only  half  of  the  corresponding  static  value.  One  element  length  away 

from  the  edge  of  the  bearing  ring,  about  5 inches,  the  peak  bearing  pressure  drops  to  i 

14,000  psi  as  Figure  5-26  shows.  A similar  though  less  drastic  drop  is  also  observed  in  | 

the  bearing  load  one  element  length  back  of  the  contact  surface.  These  results  indicate 
that  the  bearing  load  is  highly  localized.  The  fact  that  the  bearing  ring  extends  into 
and  forms  an  integral  part  of  the  base  of  the  frame  makes  interpretation  of  this  kind  even 
more  difficult 

Analysis  of  the  response  of  the  tubular  section  of  the  shelter  is  somewhat  limited 
because  of  the  early  cut-off  time.  However,  some  quantitative  design  values  can  still  be 
obtained  for  the  fore  portion  of  tubular  section,  l.e.,  the  region  behind  the  transition 
section. 

From  the  axial  stress  plots  of  Figures  5-27  for  sections  behind  the  headworks 
the  peak  mean  axial  compression  Is  estimated  to  he  9,000  psi.  Similarly,  the  peak  hoop 
stress  and  the  peak  in-plane  bending  can  be  obtained  from  Figure  5-28  for  the  tube  cross- 
section  Immediately  behind  the  headworks  and  from  Figure  5-29  for  the  cross-section  half 
tube  diameter  behind.  Note  that  whereas  the  crown  immediately  behind  the  headworks  is 
flattened  (the  outer  fibre  stress  is  higher  than  the  inner  fibre  stress  at  the  crown. 

Figure  5-28)  the  crown  half  tube  diameter  behind  the  headworks  shows  the  opposite  response 
(the  outer  fibre  stress  is  lower  than  the  inner  fibre  stress  at  the  crown.  Figure  5-29). 

The  difference  in  the  response  pattern  is  depicted  in  Figure  5-30.  It  should  be  noted  that  ' 

at  this  early  time  all  tubular  section  response  is  induced  by  loading  up  front  and  ground 
motion  induced  by  airblast  at  the  berm  has  yet  to  come  into  play. 
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Strains — as  might  be  expected  from  the  description  of  the  closure  deformation 
mode,  the  shear  strains  in  the  closure,  especially  for  points  near  the  support  are  signif- 
icant. A maximum  of  14  percent  in  the  steel  backplate  is  found  near  the  bottom  edge  and 
shear  strains  of  2 to  6 percent  are  not  uncommon.  A typical  strain  invariant  plot  (/l ' ^ 
versus  Ij^)  for  a point  in  the  backplate  is  given  in  Figure  5-31. 

Shear  strains  in  the  closure  concrete  are  in  general  of  the  order  of  1 to  3 per- 
cent, with  the  highest  strain  (8  percent)  occurring  in  the  vicinity  of  the  top  and  bottom 
shear  crack.  A typical  strain  invariant  plot  is  given  in  Figure  5-32  which  is  based  on 
the  strain/time  histories  of  Figure  5-33  for  a point  near  the  upper  support.  Comparing 
Figure  5-33  to  Figure  5-34  which  is  for  a point  at  the  center  of  the  closure  clearly 
indicates  that  much  deformation  occurs  near  the  support  while  the  center  part  of  the  clo- 
sure is  relatively  unstrained. 

5.2.  CALCULATION  2B— THREE-DIMF.NSIONAL  CALCULATION  WITH  HEAP-ON  INCIDENCE;  INELASTIC 
MODEL  FOR  CONCRETE  ONLY 

Model  and  Alrblast  Load — the  basic  model  and  airblast  loading  are  identical  to 
those  of  Calculation  2A.  The  only  changes  in  the  model  are:  the  steel  side  and  backplates 
are  assumed  elastic  and  only  one  element  is  used  across  the  plate  thickness.  This  is  done 
for  two  reasons.  Dropping  the  number  of  elements  across  the  plate  thickness  from  two  to 
one  relaxes  the  constraint  on  the  integration  time  step  somewhat  and  computation  time  per 
major  cycle  is  decreased.  Secondly,  the  change  is  made  to  investigate  what  effect  details 
In  modeling  of  the  backplate  have  on  the  overall  response  of  the  closure  and  the  main 
shelter  structure.  Using  an  elastic  material  for  steel  will  also  show  up  the  effect  of 
steel  properties  (A36  versus  A415,  for  instance)  on  closure  response.  As  a final  change, 
the  overburden  is  given  a rather  high  stiffness  (wave  speed  of  3,500  fps)  in  order  to  il- 
lustrate the  effect  of  this  parameter  on  structural  response. 

Deformation  Mode — the  deformation  pattern  of  the  headworks  at  about  15  msec,  after 
alrblast  arrival  is  Illustrated  in  Figure  5-35.  The  displacements  are  magnified  as  in 
Calculation  2A.  The  early  time  response  is  seen  to  have  the  headworks  pushed  back  and  down; 
the  latter  is  no  doubt  due  to  the  stiff  overburden  assumed.  In  Figure  5-36  are  shown  the 
displaced  profiles  of  the  backplate  vertical  center  line  at  several  action  times.  Profiles 
for  the  horizontal  center  line  Is  given  in  Figure  5-37.  The  deformation  pattern  of  the 
closure,  that  of  punching,  is  relatively  unaffected  by  the  elastic  backplate  assumption. 

A spot  comparison  of  the  deformation  in  the  concrete  closure  indicates  that  the 
closure  response  is  also  similar  to  that  in  Calculation  2A. 
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A comparison  of  Figure  5-38  which  shows  the  closure  eg  velocity  time  history 
with  its  counterpart  from  Calculation  2A  (Figure  5—11)  also  confirms  this  view.  Minor 
differences  are  observed  in  the  downward  velocity. 

Stresses — a typical  stress  invariant  plot  for  points  in  the  concrete  closure  is 
given  in  Figure  5-39  which  is  similar  to  the  results  obtained  in  Calculation  2A.  Figure 
5-40  consists  of  stress  invariant  plots  for  points  in  the  backplate.  Because  of  the 
elastic  assumption  there  is  no  yielding.  However,  two  trends  are  evident.  Had  the  steel 
been  replaced  by  a A36  material  with  octahedral  shear  stress  limit  of  20,800  psi , inelastic 
deformation  would  certainly  occur  near  the  edges,  an  observation  which  is  confirmed  in 
Calculation  2A.  Secondly,  for  the  plate  to  remain  elastic  a steel  with  grades  much  higher 
than  A36  is  necessary.  However,  this  upgrading  does  not  seem  to  have  significant  effect 
on  the  behavior  of  the  closure  as  was  pointed  out  earlier. 

Behavior  of  the  two  corners  of  the  frame  is  also  examined.  While  the  concrete 
response  is  similar  to  that  of  Calculation  2A,  as  is  exemplified  in  Figures  5-41  and  5-42, 
the  elastic  assumption  for  steel  brings  out  an  interesting  observation.  Referring  to 
Figures  5-43  and  5-44  for  points  in  the  upper  and  lower  steel  frame  corners,  respectively, 
one  sees  now  that  the  lower  corner  sustains  a higher  shear  stress  whereas  the  upper  corner 
is  lightly  stressed  by  comparison.  Limiting  the  shear  load  in  the  steel  transfers  part 
of  the  bottom  load  to  the  top  as  observed  in  Calculation  2A  and  pushes  both  the  upper 
and  lower  corner  beyond  the  elastic  range. 

From  invariant  plots  obtained  for  points  in  the  headworks  and  tubular  section,  it 
is  observed  that  the  response  is  quite  similar  to  that  of  Calculation  2A. 

Loads  exerted  by  the  closure  on  the  bearing  ring  are  shown  in  Figure  5-45.  Since 
the  direct  hearing  (contact)  area  is  modeled  by  one  element  width  the  contact  stress  given 
in  the  figure  lies  roughly  between  the  upper  and  lower  bound  values  given  by  Figures  5-25 
and  52-6  for  Calculation  2A.  Again  the  bearing  load  transfer  does  not  seem  to  be  signifi- 
cantly effected  by  the  lack  of  details  in  modeling  of  the  backplate  in  this  case. 

Similar  remarks  can  be  made  about  the  propagation  of  the  longitudinal  stress  along 
the  shelter.  One  significant  difference  is  noted,  however;  referring  to  Figure  5-46 
which  shows  the  longitudinal  stress/time  histories  for  points  in  the  tubular  sections 
behind  the  headworks  the  differential  between  the  longitudinal  stresses  at  the  crown 
and  Invert  (a  measure  of  the  bending  moment  acting  in  the  tubular  section)  is  larger 
than  that  shown  in  Figure  5-27  f(  r Calculation  2A;  furthermore,  it  does  not  seem  to 
diminish  with  time  or  distance  along  the  tube.  This  bending  moment  can  be  attributed  to 
the  effect  of  the  downslap  force  acting  on  the  headworks.  The  stiff  backfill  properties 
assumed  make  the  effect  far  more  significant  than  it  is  In  Calculation  2A. 


The  in-plane  bending  moment  indicates  an  ovaling  deformation  of  the  circular 
cross-section  immediately  behind  the  transition  oection  and  of  higher  mode  deformation 
for  the  section  about  6 feet  back  (see  Figure  5-47),  not  unlike  the  results  obtained  in 
Calculation  2A. 

5.3.  CALCULATION  2C~THREE-DIMENSI0NAL  ELASTIC  CALCULATION  WITH  HEAD-ON  INCIDENCE 

Model  and  Airblast  Load — the  model  and  airblast  load  of  Calculation  2B  are  used 
with  all  material  properties  assumed  elastic.  Consequently,  the  response  of  the  closure 
so  obtained  corresponds  to  one  which  is  designed  to  remain  elastic  under  the  prescribed 
load. 

Deformation  Mode — the  deformation  patterns  of  the  headworks  at  15  msec,  and 
22  msec. , after  airblast  arrival  are  Illustrated  in  Figure  5-48,  the  displacements  having 
been  magnified.  Aside  from  the  expected  backward  and  downward  movement,  there  is  an  al- 
most rigid-body  like  appearance  about  the  deformation  pattern  which  is  typical  of  an 
elastic  solution  (see  Phase  1 results,  for  instance).  Note  that  the  closure  support  point 
is  no  longer  pushed  straight  backward  as  is  observed  in  Calculations  2A  and  2B.  Instead 
the  support  point  moves  as  part  of  the  rigid  body  which  is  the  headworks. 

In  Figure  5-49  are  shown  the  displaced  profiles  of  the  backplate  vertical  center 
line  at  several  action  times.  Profiles  for  the  horizontal  center  line  are  given  in 
Figure  5-50.  Maximum  deflection  at  mid-span  relative  to  the  support  is  about  two  inches. 

The  velocity-time  history  of  the  closure  eg  is  shown  in  Figure  5-51  and  is  os- 
cillatory as  typical  of  an  elastic  system.  The  dominant  period  of  the  longitudinal  com- 
ponent is  between  7 and  8 msec.,  giving  the  closure  a natural  frequency  of  125  to  140  Hz. 
The  downward  velocity  on  the  other  hand  contains  a 250  Hz  component  superimposed  on  the 
overall  headworks  motion. 

Stresses — typical  stress  invariant  plots  for  points  in  the  concrete  closure  are 
given  in  Figure  5-52.  The  concrete  near  the  support  edges  would  have  yielded  if  the  yield 
criterion  as  shown  also  on  the  figures  had  been  imposed.  The  concrete  at  the  center  of  the 
closure,  on  the  other  hand,  remains  elastic;  it  would  appear  to  remain  so  even  when  the 
yield  criterion  is  Imposed. 

The  behavior  of  the  backplate  is  illustrated  in  the  Invariant  plots  of  Figure  5-53 
for  points  near  the  support  edges  and  for  a point  at  the  center  of  the  backplate.  As 
expected  the  plate  center  sustains  high  tens*le  stresses  as  the  closure  deforms  as  a plate 
under  uniform  load. 
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Behavior  at  the  upper  and  lower  corners  of  the  frame  is  shown  In  Figures  5-54 
through  5-57.  It  Is  seen  that  both  steel  and  concrete  (being  elastic)  sustain  tensile 
stresses.  One  can  readily  surmise  that  if  the  shear  and  tensile  stress  bearing  capability 
of  the  concrete  is  limited  most  of  the  tensile  stress  will  be  transferred  to  the  steel 
liner  and  this  is  indeed  the  case  as  shown  in  Calculation  2B  (Figures  5-43  and  5-44), 

Stress  invariant  plots  for  a point  in  the  tubular  section  of  the  shelter  are 
shown  in  Figures  5-58  and  5-59.  When  compared  with  the  yield  criterion  for  concrete  they 
show  that  inelastic  deformation  will  take  place  at  the  crown  and  springline  and  to  a 
nominal  extent  at  the  invert. 

Loads  exerted  by  the  closure  on  the  bearing  ring  are  shown  in  Figure  5-60.  Note 
that  the  top  and  side  bearing  pressures  are  almost  identical  and  are  higher  than  that  at 
the  bottom,  a behavior  also  observed  in  the  elastic  results  of  Phase  1.  Further  the  time 
histories  have  a dominant  oscillatory  component  with  a frequency  close  to  that  of  the  clo- 
sure vibration.  The  peak  bearing  load  is  about  32,000  psl. 

Longitudinal  transfer  from  headworks  into  the  tubular  section  follows  the  usual 
pattern  with  an  amplification  factor  of  slightly  higher  than  two  (see  Figure  5-61).  The 
axial  stresses  at  the  crown  and  springline  are  higher  than  that  at  the  invert.  It  is  also 
worthwhile  to  point  out  that  the  pulse  behaves  like  an  exponential  pulse  with  a rise  time 
of  5 msec,  (the  alrblast  load  rise  time)  and  main  pulse  duration  of  20  msec. 
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Figure  5-Aa.  Soil  model  in  uniaxial  strain,  axial  stress  versus  axial  strain 


Soil  model  in  uniaxial  strain,  stress  difference  versus  pressure 
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Figure  5-4c.  Comparison  of  backfill  models,  uniaxial  stress/strain  curves 


Figure  5-5.  Comparison  of  overpressure  waveforms  at  front  face,  front-on  incidence.  Calculations  2A 
2B  and  2C 


f'igure  5-7.  Deformation  mode  of  backplate,  vertical  centerline.  Calculation  2A 


ri^ontal  centerline.  Calculation  2A 


Inelastic  Zone 


5-9a.  Development  of  inelastic  zones  in  backplate.  Calculation  2A.  Arlval  time  of  airblast  is 
113  msec. 


5-9b.  Development  of  Inelas 


Inelastic  Zone 


airblast  is  113  msec. 


Oc.  Development  of  Inelastic  zone,  in  closure  concrete.  Calculation  2A  (continued).  Arrival  time 
of  airblast  is  113  msec. 


Figure  5-lOd.  Development  of  inelastic  zone,  in  closure  concrete.  Calculation  2A  (continued).  Arrival 
time  of  airblast  is  113  msec. 


-lOe.  Development  of  inelastic  zone,  in  closure  concrete,  Calculation  2A  (continued).  Arrival 
time  of  alrblast  is  113  msec. 


figure  5-12.  Typical  stress  invariant  plot  for  points  In  closure  concrete  near  support.  Calculation  2A 


points  In  center  of  closure  concrete.  Calculation  2A 


Figure  5-14.  Typical  stress  Invariant  plot  for  points  In  backplate  near  support.  Calculation  2A 


unsupported,  Calcula 


Invariant  plot  for  upper  frame  corner,  outer  steel  lining.  Calculation  2A 


ress  Invariant  plot  for  lower  frame  corner,  inner  steel  lining.  Calculation  2A 


Stress  invariant  plot  for  lower  frame  corner,  outer  steel  lining. 


Figure  5-18,  Typical  stress  invariant  plot  for  concrete  in  upper  frame  corner,  Calculation  2A 


Typical  stress  invariant  plot  for  concrete  in  lower  frame  corner.  Calculation  2A 


Figure  5-20.  Typical  stress  Invariant  plot  for  concrete  In  frame  gussets.  Calculation  2A 
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Figure  5-28.  Hoop  stress/time  histories  for  tube  section  behind  headworks.  Calculation  2A 


invariant  plot  for  points  in  backplate  near  support.  Calculation  2A 


1-33.  Typical  strain/timc  histories  for  concrete  in  closure,  near  support,  Calcu 
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Figure  5-34.  Typical  strain/tlme  histories  for  concrete  in  closure,  near  center.  Calculation  2A 


Figure  5-36.  Deform.Tt Ion  modes  of  backplate  vertical  centerline,  Calculation  2B 
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Figure  5-46.  Longitudinal  stress/t ime  histories,  tube  sections  behind  headworks.  Calculation  2B 
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Figure  5-52a.  Typical  stress  Invariant-plot  for  points  In  closure  concrete,  near  support,  Calculation  2C 


Figure  5-54.  Typical  stress  invariant  plot  for  concrete  in  upper  frame  corner.  Calculation  2C 


Typic 


Typical  stress  Invariant  plot  for  concrete  in  tubular  section,  crown  and  sprlngllne 
Calculation  2C 


Figure  5-59.  Typical  stress  invariant  plot  for  concrete  in  tubular  section,  invert.  Calculation  2C 


histori 


Figure  5-61b. 


MATERIAL  PARAMETERS  FOR  S4  MODEL 


*Used  in  conjunction  with  different  tension  cutoffs  (7.2  ksi  to  0.5  ksi)  to  model  various  reenforced 
concrete  in  the  structure,  headwall,  slab  floor  and  grout. 


SECTION  VI 


INTERPRETATION  OF  PHASE  2 RESULTS 

Major  results  of  Calculations  2A,  2B  and  2C  all  based  on  the  S4  shelter  conflgura- 
; tion  and  head-on  incidence  are  discussed  below.  Since  emphasis  of  the  study  is  not  on  struc- 

■ tural  response  per  se,  the  discussion  on  structural  response  will  be  kept  to  a minimum;  de- 

tails are  given  only  to  provide  a common  data  base  for  subsequent  discussions  which  address 
the  main  purpose  of  the  study. 

j 6.1.  DEFORMATION  MODE 

1 

t The  overall  deformation  pattern  of  the  shelter  is  consistent  with  that  observed 

in  Phase  1 which  is  based  on  the  SI  configuration.  The  dominant  force  comes  from  the  direct 
head-on  load  acting  on  the  front  face.  The  structure  responds  like  a beam-column  supported 
by  the  soil  medium;  the  headworks  is  pushed  back  in  the  longitudinal  direction  by  tliis  load 
with  some  downward  motion  due  to  the  inclined  front-face  and  overburden  load.  The  motion  of 
the  headworks  in  turn  induces  some  bending  of  the  tubular  section  behind  the  headworks,  not- 
ably in  the  longitudinal  plane.  In-plane  bending  of  the  tubular  section  of  the  shelter  also 
exists  as  a result  of  the  longitudinal  load  at  first  and  later  on  as  a result  of  the  over- 
burden load  acting  on  that  section. 

An  imbalance  of  the  longitudinal  stress  exists  in  the  upper  (top)  and  lower  (bottom) 
portions  of  the  headworks.  The  imbalance  is  due  to  two  factors;  (a)  the  shielding  effect  of 
the  soil  which  protects  the  lower  portion  of  the  front  face  from  direct  airblast  load;  and 
(b)  the  higher  load  bearing  capability  of  the  medium  supporting  the  lower  portion  of  the 
shelter.  Consequently,  the  top  portion  of  the  headworks  and  part  of  the  tubular  section 
which  follows  sustain  greater  motion  than  the  lower  portion. 

The  downward  motion  of  the  headworks  on  the  other  hand  depends  on  the  overburden 
load  and  the  timing  of  this  downslap  depends  on  the  wave  propagation  speed  of  the  berm 
medium.  For  a soft  berm  (loose  backfill)  the  overburden  load  does  not  reach  the  structure 
for  the  first  10  msec,  or  so  after  airblast  arrival  at  the  front  face;  it  is  effectively 
screened  out  for  short-time  response  consideration.  Little  downward  motion  of  the  headworks 
results.  When  the  berm  material  has  a higher  wave  speed  (dense  backfill),  the  downslap 
f rce  Is  Immediate  and  downward  motion  is  almost  concurrent  with  the  longitudinal  motion. 

Slme  the  bulk  of  the  headworks  remains  elastic  it  is  not  surprising  that  the 
•i.f,;  Jetorm.it  i>'n  pattern  as  summarized  above  holds  true  for  Calculations  2A,  2B  and  2C. 

: I,  - !•  deform.it  Ions  in  the  frame  and  closure  support  areas  account  for  some  of  the 
I.  t«i!i  In  this  general  pattern  (see  Figures  5-35  and  5-48  for  instance). 
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The  closure  Is  pushed  straight  back  by  the  incident  airblast  into  the  bearing  ring 
and  its  support  area  in  the  headworks.  Consequently,  the  steel  bearing  ring  and  in  particular 
the  concrete  region  behind  it  are  severely  deformed.  The  10*^  incline  of  the  closure  also 
causes  a downward  component  of  the  closure  motion  which  is  evident  at  the  bottom  edge  of  the 
closure  and  lower  bearing  area. 

The  deformation  of  the  closure  is  very  much  affected  by  the  closure  material  des- 
cription. This  is  because  the  closure  behaves  like  a thick  plate  in  flexure  (which  is  slight 
because  of  the  large  thlckness-to-span  ratio)  in  the  elastic  range  and  like  a slab  in  punch- 
ing deformation  in  the  inelastic  range.  The  two  deformation  modes  are  significantly  differ- 
ent as  can  be  seen  by  comparing  Figures  5-7  and  5-49  showing  the  deformed  profile  of  the 
backplate.  Because  of  its  importance  in  subsequent  discussion  the  extrusion  or  punching  mode 
will  be  described  in  the  following  paragraphs. 

For  the  closure  the  airblast  load  at  the  front  face  (lateral  surface  of  the  slab) 
is  transformed  into  shear  forces  in  the  slab.  High  shear  stresses  first  develop  in  the  con- 
crete portion  near  the  support  edge  as  a result  of  the  high  frontal  load  and  the  equally 
high  reaction  forces  along  the  support  edge.  This  shear  stress  is  highest  at  the  bottom 
concrete  which  is  next  to  the  backplate.  Upon  yielding  of  this  bottom  concrete,  i.e., 
when  the  shear  bearing  capability  of  the  concrete  is  reached,  the  surplus  shear  stress  is 
transferred  to  the  concrete  immediately  in  front  which  in  turn  yields.  This  process  which 
continues  for  sufficiently  high  applied  load  until  the  front-most  concrete  in  the  closure 
is  similarly  affected,  is  often  referred  to  as  the  shear  crack  propagation  phenomenon. 

Once  the  propagation  is  complete,  i.e.,  it  extends  through  the  thickness  of  the 
closure  concrete  any  additional  load  must  be  carried  by  the  steel  backplate.  When  the 
additional  load  is  sufficiently  high  as  is  the  case  in  Calculations  2A  and  2B  the  backplate 
undergoes  a shear  deformation  at  the  support  edges.  The  closure  in  effect  is  extruded  or 
punched  into  the  headworks  cavity  shearing  at  the  edges,  but  with  the  overall  backplate  re- 
maining relatively  flat. 

This  punching  deformation  of  the  closure  is  not  sensitive  to  the  backplate  details 
at  least  as  far  as  the  first  10  to  15  msec,  of  response  is  concerned.  Changing  the  backplate 
material  from  elastic  to  elasto-plastic  and  increasing  the  number  of  elements  across  the 
plate  thickness  from  I to  2 do  not  seem  to  have  significant  effect  on  the  deformation  pattern 
of  the  closure.  What  is  found  to  have  significant  influence  on  the  closure  response  is  the 
shear  bearing  capability  of  the  closure  concrete.  Changing  the  closure  concrete  material 
from  elastic  (Infinite  shear  bearing  capability)  to  elasto-plastic  (limited  shear  bearing 
capability  according  to  confinement  pressure)  causes  an  estimated  four-fold  increase  In  the 
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maximum  deflection  as  measured  by  the  center  deflection  of  the  backplate  relative  to  its 
support. 

This  sensitivity  can  also  be  observed  by  examining  the  motion  of  the  closure  eg 
for  the  two  cases,  elastic  closure  and  inelastic  closure,  reproduced  in  Figures  6-la  and 

6-lb.  Discounting  the  10°  incline,  the  Y-component  of  the  velocity  can  be  viewed  as  the 

component  perpendicular  to  the  plane  of  the  closure.  In  the  elastic  case  the  eg  is  seen  to 
exhibit  an  oscillatory  response  with  a period  of  about  8 msec.  The  velocity  is  similar  to 
that  of  a single  degree-cf-f reedom  elastic  system  subjected  to  a triangular  pulse  with  dura- 
tion several  times  that  of  the  period  of  the  system.  The  onset  of  inelasticity  and  the 

formation  of  shear  cracks,  however,  change  all  that.  The  velocity-time  history  of  the  in- 

elastic closure  eg  is  not  unlike  that  of  a single  degree-of-f reedom  elasto-plastic  system. 

These  conclusions  have  great  ramifications  in  cost-hardness  trade-off  of  the 
shelter  (closure).  High  cost  and  low  sensitivity  for  the  steel  backplate  suggest  that  it 
may  not  merit  the  amount  of  attention  that  it  has  received.  On  the  other  hand,  low  cost 
plus  great  sensitivity  for  the  closure  concrete  dictates  that  it  be  considered  as  a very 
important  cost-hardness  trade-off  parameter. 

Once  the  closure  is  extended  Into  Its  punching  deformation  mode  the  airblast 
impulse  in  the  main  pulse  duration  governs  the  response,  i.e.,  how  fast  or  slow  the  pres- 
sure pulse  decays  from  its  peak  value.  Hence,  it  is  expected  that  higher  frequency  contents 
which  may  be  important  to  the  closure  response  in  the  elastic  realm  (for  vibration  and 
resonance  considerations)  are  no  longer  as  important  since  they  do  not  greatly  change  the 
load  impulse.  In  other  words  if  the  closure  is  designed  to  respond  in  the  punching  mode 
the  only  overpressure  wave  form  parameter  of  Interest  (aside  from  peak  magnitude)  is  the 
total  Impulse  and  how  fast  this  Impulse  is  achieved. 

The  tubular  section  of  the  shelter  immediately  behind  the  transition  section 
undergoes  an  ovaling  deformation  partly  due  to  the  oncrushing  headworks  (short-term  effect) 
and  partly  due  to  the  overburden  load  from  the  top  (longer-term  phenomenon).  The  section 
about  half  shelter  diameter  behind  the  headworks  shows  the  deformation  pattern  of  a mode 
higher  than  oval Ing, poss lb ly  the  n = 3 mode  although  the  results  obtained  are  not  compre- 
hensive enough  to  be  conclusive. 

6.2.  STRESS  DISTRIBUTION  IN  SHELTER 

The  airblast  load  acting  on  the  front  face  (frame/closure)  is  transmitted 
through  the  headworks  to  the  tubular  section.  The  peak  magnitude  of  this  longitudinal  stress 
is  amplified  in  going  from  the  headworks  to  the  tubular  section  because  of  the  abrupt  drop 
in  load  bearing  area  available  in  the  latter  although  the  amount  of  amplification  is  not  as 
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high  as  would  be  expected  based  on  bearing  area  consideration  alone.  There  are  two  contribut- 
ing factors  for  this:  the  load  bearing  limit  in  the  tubular  section  which  corresponds  to 
the  strength  of  concrete  used  (6,000  psi,  unconfined),  and  the  load  transfer  between  closure 
and  bearing  plate,  which  will  be  discussed  presently. 

Accompanying  the  amplification  in  longitudinal  stress  is  a change  in  wave  form 
partly  due  to  geometric  dispersion  and  partly  due  to  Inelastic  deformation.  The  change 
due  to  dispersion  can  be  measured  by  the  elastic  results  of  Calculation  2C  and  the  contri- 
bution due  to  inelastic  deformation  by  comparing  these  v ith  the  results  of  Calculation  2A 
or  Calculation  2B. 

The  longitudinal  stresses  in  the  upper  portion  of  the  shelter  are  invariably 
higher  than  those  in  the  lower  portion  as  has  been  mentioned  in  discussing  the  deformation 
mode.  This  difference  in  the  tubular  section  tends  to  diminish  farther  away  from  the  head- 
works.  Although  not  verified  directly  it  is  expected  that  the  longitudinal  stress  differ- 
ential vanishes  altogether  about  several  tube  diameters  down  the  tube. 

The  axial  stress  transmitted  to  the  tubular  section  is  sufficiently  high  in  areas 
close  to  the  headworks  that  inelastic  deformation  is  observed  there  both  at  the  crown  and 
springlines  and  to  a lesser  extent  at  the  Invert. 

The  stress  state  in  the  closure  itself  has  also  been  covered  in  connection  with 
the  discussion  on  the  punching  deformation  of  the  closure.  In  the  following  paragraphs 
the  stress  state  at  the  closure/bearing  Interface  will  be  described. 

Much  like  a plate  under  lateral  load,  the  airblast  loading  acting  on  the  closure 
induces  reaction  forces  at  the  bearing  ring.  In  tlie  static  situation  the  ratio  of  the  re- 
action pressure  to  the  uniform  applied  load  is  proportional  to  the  ratio  of  the  loaded  sur- 
face to  the  supported  area.  Results  of  Phase  2 indicate  that  the  dynamic  reaction  load 
Induced  in  the  bearing  ring  is  lower  than  the  static  equivalent.  Furthermore,  the  reaction 
load  is  concentrated  at  the  innermost  edge  of  the  ring  which  is  consistent  with  the  punch- 
ing deformation  of  the  closure.  The  average  bearing  pressure  is  found  to  be  insensitive  to 
the  details  of  the  Interface  model  whether  the  backplate  is  elastic  or  inelastic  and  whether 
it  has  two  elements  or  only  one  element  across  the  plate  thickness.  A fine  mesh  representa- 
tion of  course  should  be  used  to  represent  the  bearing  area  if  the  stress  concentration 
nature  of  the  bearing  load  is  to  be  reproduced. 

The  model  used  in  Phase  2 assumes  a bonded  Interface  between  the  closure  and  bear- 
ing ring.  The  actual  physical  interface  condition  is  of  course  complicated  and  appears 
time  dependent.  The  gap  pressure  wave  which  precedes  the  main  airblast  load  at  the  front 
face  will  tend  to  slip  past  the  bearing  Interface  into  the  headworks  cavity.  This  was 
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observed  In  the  HAVE  HOST  SI  test  (Reference  11)  and  the  so-called  blow-by  is  sometimes  so  > 

significant  that  it  knocks  off  equipment  and  cables  inside  the  headworks.  During  this  period  i 

which  occurs  when  the  closure  is  relatively  unloaded  the  interface  condition  is  one  of  slip.  ; 

Upon  the  establishment  of  the  full  airblast  load  on  the  closure  the  interface  is  one  of  ~ 

friction  or  no-slip  depending  on  the  normal  force  and  the  relative  motion  of  the  backplate 
with  respect  to  the  bearing  ring.  When  loaded  on  its  front  face  the  closure  will  tend  to  • 

extend  sideways.  The  confinement  pressure  on  the  other  hand  will  tend  to  move  it  in  the 
opposite  direction.  What  actually  takes  place  is  not  totally  clear  and  no  doubt  depends  on 
the  relative  magnitude  of  the  confinement  (gap)  pressure  and  front  airblast  load.  The  effect 
of  interface  condition  on  shelter  response  is  expected  to  be  minor.  The  interface  (or  edge) 
condition  may  be  important  to  the  response  of  th*.  closure , and  its  effect  needs  to  be  assessed. 


The  aformentioned  aspects  of  stress  distribution  in  the  shelter,  namely,  the  bear- 
ing pressure,  the  peak  mean  longitudinal  stress  and  bending  moment  in  the  tubular  section, 
etc. , are  found  to  be  insensitive  to  the  modeling  details  of  the  closure  (comparing  results 
of  Calculations  2A  and  2B) . Although  the  observation  is  based  on  short-term  response  we 
believe  that  it  applies  also  to  long-term  response  as  well,  leading  one  to  conclude  that  as 
far  as  the  overall  response  of  the  shelter  is  concerned  a crude  model  of  the  closure  suf- 
fices. By  the  same  token  for  detailed  analyses  of  the  closure  the  headwork  and  tubular 
section  may  be  replaced  by  some  equivalent,  but  much  simplified  system  suitable  chosen  to 
maintain  the  correct  impedance  and  energy  transmission  properties. 

Based  on  results  obtained  in  Phase  2 and  the  discussions  of  Sections  6.1  and 
6.2,  the  following  observations  on  the  effect  of  inelasticity  (and  elasticity)  on  shelter 
response  can  be  made. 


A.  Effect  of  Inelasticity 

The  foremost  effect  of  inelasticity  on  shelter  response  is  manifested  in  the 
deformation  and  response  of  the  closure.  Limiting  the  shear  bearing  capability  of  the 
closure  concrete  results  in  the  formation  of  a shear  crack  causing  the  closure  to  respond 
in  the  punching  mode  wherein  it  is  extruded  into  the  headworks  cavity  for  sufficiently 
high  loadings.  Consequently,  the  steel  bearing  ring  and  the  concrete  region  behind  it  are 
also  severely  stressed.  This  distress  is  more  evident  at  the  bottom  edge  of  the  closure 
and  lower  bearing  area. 

An  indirect  effect  of  Inelasticity  in  the  closure  is  that,  once  the  closure  is 
in  the  punching  mode,  the  closure  deformation  is  governed  by  the  loading  impulse  and  is  no 
longer  as  sensitive  to  high  frequency  components  of  the  overpressure  wave  form  as  is  when 
it  remains  elastic. 
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Inelastic  deformations  are  found  in  the  frame,  notably  at  the  upper  and  lower 
inner  corners.  The  effect  of  inelasticity  is  local  in  the  sense  that  it  does  not  affect 
the  response  of  the  main  headworks  and  shelter  structure. 

Incipient  plastic  deformations  are  also  found  in  the  tubular  section  immediately 
behind  the  headworks  where  the  high  axial  load  at  the  crown  is  limited  by  the  concrete 
strength  and  is  redistributed  to  the  lower  portion  of  the  tube  causing  a more  uniform  dis- 
tribution over  the  tube  cross-section. 

B . Elastic  Versus  Inelastic  Analysis 

By  comparing  the  elastic  and  inelastic  response  obtained  it  is  observed  that  an 
elastic  analysis  of  the  shelter  is  valuable  in  identifying  both  the  overall  deformation 
modes  of  the  structure  especially  those  that  are  uneffected  by  the  onset  of  inelasticity 
and  possible  regions  of  distress.  Some  examples  from  the  S4  configuration  are  the  head- 
works  deformation,  load  transfer  from  closure  to  headworks,  load  transfer  from  headworks 
to  tubular  section,  stresses  in  the  frame,  stresses  in  the  closure  and  tubular  section, 
and  response  of  the  closure  provided  that  it  remains  elastic. 

The  usefulness  of  an  elastic  analysis  can  be  further  enhanced  when  its  results 
are  interpreted  based  on  sound  engineering  Judgment.  For  example,  by  comparing  the 
elastic  stresses  against  the  failure  criterion  of  the  material,  an  indication  of  whether 
the  material  may  be  overstressed  can  be  obtained.  Imposing  a limit  on  the  load  bearing 
capability  of  the  material  an  experienced  analyst  might  he  able  to  transfer  the  surplus 
load  to  other  parts  of  the  structure  as  is  indeed  what  actually  happens.  The  transfer  be- 
tween the  upper  and  lower  corner  of  the  frame,  the  transfer  between  the  upper  and  lower 
portion  of  the  tubular  sections,  the  transfer  from  rear  closure  concrete  to  the 
front  are  some  of  the  lessons  learned  in  the  case  studies  of  Phase  2. 

Elastic  analyses  of  the  shelter  can  provide  useful  transient  (short  term) 
response  data  here  because  the  structure  under  study  is  semi-buried  and  subjected  to  a 
direct  loading  which  dominates  the  short-term  response.  The  elastic  analyses  will  be 
more  effective  if  some  allowance  is  made  for  energy  dissipation  inherent  in  the  real 
material . 

For  a detailed  quantitative  analysis  of  the  shelter  and  particularly  the  clo- 
sure, there  is  no  substitute  for  an  inelastic  analysis.  Whereas  the  motion  of  the  clo- 
sure eg  can  be  approximated  readily  by  an  equivalent  single  degree-of-f reedom  system, 
the  punching  deformation  mode  and  the  shear  crack  propagation  phenomenon  are  prevalent 
only  in  an  inelastic  three-dimensional  dynamic  analysis  of  the  structural  component 
involved. 
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Figure  6-1.  Comparison  of  velocity/time  histories  of  closure  eg 


SECTION  VII 


CONCLUSIONS 

This  section  begins  with  a summary  of  all  modeling  considerations  found  to  have 
significant  impact  on  shelter  response.  Basic  deformation  modes  of  the  shelter  under  dif- 
ferent loading  conditions  are  described  and  their  sensitivity  to  airblast  loading  charac- 
teristics summarized.  The  effect  of  inelasticity  on  modeling,  deformation  pattern  and  air- 
blast  sensitivity  is  then  delineated.  Once  the  effect  of  material  nonlinearities  has  been 
identified  the  merit  and  limit  of  the  use  of  elastic  analyses  (and  hence  inelastic  analyses) 
will  become  apparent. 

With  knowledge  gleaned  from  these  conclusions  a viable  approach  to  performing 
dynamic  analyses  of  the  shelter  is  developed  and  described.  The  procedure  involves  model- 
ing the  shelter-like  structure  and  its  medium  according  to  the  guidelines  established, 
factoring  in  those  characteristics  of  the  air  overpressure  which  have  been  found  signif- 
icant and  separating  the  analysis  into  complementary  components;  short-term  analysis 
and  long-term  analysis. 

7.1.  STRUCTURAL  MODELING  CONSIDERATIONS 

Several  regions  of  the  structure  need  to  be  modeled  with  attention.  These 
are  the  closure,  the  front  face,  the  inside  corners  in  the  frame,  the  gap  between  the 
closure  and  frame  and  the  transition  region  including  the  front  portion  of  the  tube. 

7.1.1.  Closure 

Details  in  the  closure  model  do  not  seem  to  significantly  effect  the  response 
of  the  remainder  of  the  structure.  When  only  the  gross  response  of  the  structure  is  de- 
sired a crude  model  may  be  used  for  the  closure  in  order  to  accommodate  a large  integra- 
tion time  step  for  an  explicit  code.  This  is  essential  in  computing  the  long  time  response 
of  the  shelter.  When  accurate  response  of  the  closure  itself  is  desired  the  closure  model 
should  be  sufficiently  refined  to  preserve  its  response  characteristics  (natural  vibration 
in  the  clastic  realm  and  punching  or  shear  failure  in  the  inelastic  range).  The  remainder 
of  the  shelter,  however,  may  be  replaced  by  an  equivalent  system  provided  the  replacement 
system  retains  the  same  Impedance  and  radiation  damping. 

Closure  resporse,  in  addition.  Is  very  sensitive  to  the  material  description 
of  the  closure  constituents  and  in  the  case  of  a composite  slab,  the  sensitivity  to  con- 
crete properties  is  much  higher  than  that  to  steel.  Hence,  an  accurate  material  model 
for  the  closure  concrete  should  be  used. 
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The  effect  of  support  condition  on  closure  response  needs  to  be  assessed  but  does 
not  appear  significant.  Hence  the  closure/bearing  interface  model  may  be  crude  as  far  as 
closure  response  or  headworks  response  is  concerned.  If  it  should  become  desirable  to  ob- 
tain the  peak  bearing  load,  however,  a fine  mesh  is  necessary  since  this  load  is  concentrat- 
ed near  the  edges  and  decays  quite  rapidly  with  distance  away  from  the  edge. 

7.1.2.  Frame  and  Gap 

The  effect  of  the  gap  pressure  is  significant  especially  when  the  inner  frame  is 
unprotected.  The  pressure  tends  to  push  the  frame  outward  and  this  action  is  most  severe 
at  the  interior  corners.  Accordingly  the  gap  and  the  frame  corner  areas  should  be  included 
in  a detailed  model.  The  effect  is,  however,  local  in  that  the  headworks  response  is  un- 
affected by  it. 

The  degree  of  refinement  in  the  front  face  (frame/gap  and  closure)  should  be 
compatible  with  the  overpressure  distribution  since  the  latter  dominates  the  shelter 
response  and  bending  of  the  headworks  in  particular. 

7.1.3.  Headworks 

The  sensitive  region  in  the  headworks  appears  to  be  the  region  behind  the  bearing 
ring  which  is  severely  stressed  as  the  closure  is  pushed  against  the  support.  Otherwise 
the  bulk  of  the  headworks  is  relatively  unstressed  and  a nominal  model  suffices. 

7.1.4.  Transition  and  Tube 

A large  axial  stress  gradiant  exists  in  the  transition  region  from  headworks  to 
tube.  This  together  with  the  abrupt  change  in  cross-sectional  shape  (at  least  in  the  con- 
figurations studied)  make  the  transition  region  one  of  the  most  Important  segments  of  the 
shelter  as  far  as  overall  shelter  response  is  concerned. 

The  front  tubular  section  of  the  shelter,  i.e.,  the  portion  of  the  tube  immedi- 
ately behind  the  transition  region,  is  another  sensitive  region  in  the  shelter  response. 

In  addition  to  transmission  and  merging  of  the  axial  loads  from  the  headworks  it  transfers 
the  bending  moment  exerted  by  the  headworks  on  the  tube.  Its  deformation  is  further  con- 
strained by  the  relatively  rigid  transition  region  and  includes  contributions  from  modes 
higher  than  simple  ovaling  and  breathing  which  are  standard  for  the  remainder  of  the  tube. 

7.2.  DEFORMATION  MODE 

The  deformation  mode  under  head-on  loading  is  studied  quite  extensively  using 
elastic  and  inelastic  calculations.  The  studies  of  slde-on  and  oblique  front-on  loadings 
are  limited  and  based  exclusively  on  elastic  analyses.  The  following  discussions  are 
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limited  to  the  short-term  response  for  the  first  10  to  15  msec,  after  airblast  first 
arrival. 

7.2.1.  Head-On  Loading 

The  basic  deformation  of  the  structure  under  head-on  loading  is  one  of  axial 
compression  with  minor  variations  in  the  form  of  headworks  rotation,  longitudinal  bending 
and  downslap  motion  depending  on  details  of  the  air  overpressure  distribution  on  the  expos- 
ed surfaces,  shelter  front-face  configuration  and  berm  material. 

The  response  is  dominated  by  the  tremendous  front  load  acting  at  the  frame/ 
closure;  consequently,  is  influenced  by  its  distribution  at  the  front  face.  More  will 
be  said  about  this  in  the  section  on  airblast  characteristics.  Suffice  it  to  say  at  this 
point  that  the  headworks  is  pushed  back  into  the  tube  section  causing  severe  distress  in 
the  transition  section  and  the  front  portion  of  the  tube.  In  addition  there  is  a down- 
ward motion  induced  by  the  10°  Incline  of  the  front  face  such  as  in  the  S4  configuration 
and/or  the  overburden  load.  There  is  also  a rotation/bending  of  the  headworks,  the  nature 
and  severity  of  which  depend  on  several  factors:  any  imbalance  in  the  front  loading,  off- 
set between  center  of  pressure  and  center  of  structure  stiffness,  magnitude  of  overburden 
load  and  inclination  of  the  front  face. 

In  the  environment  considered  the  closure  deforms  in  the  punching  mode  where  it 
is  pushed  back  against  the  bearing  ring  and  extruded  into  the  headworks  cavity.  The  re- 
mainder of  the  slab  remains  relatively  flat.  Large  shear  deformations  occur  in  the  clo- 
sure both  in  the  concrete  and  steel  backplate,  mostly  along  the  support  edges.  The  shear 
deformation  initiates  at  the  back  concrete,  then  propagates  toward  the  front  face.  Upon 
completion  of  the  "shear  cone,"  the  excess  load  is  transferred  to  the  steel  backplate 
which  then  yields  in  the  punching  mode. 

Cap  pressure  which  acts  in  the  opening  between  the  closure  and  frame  tends  to 
push  out  on  the  frame  and  contain  the  closure.  As  a consequence  the  frame  sustains 
high  hoop  tension  especially  at  the  two  Inner  corners  where  tension  failure  at  the  con- 
crete and  yielding  of  the  steel  lining  is  observed.  The  gap  is  also  opened  by  this  in- 
ternal pressure  and  the  lower  lip  opens  more  than  the  upper  lip. 

The  high  axial  stress  and  any  longitudinal  bending  of  the  headworks  are  trans- 
mitted to  the  tubular  section  with  the  fore  portion  of  the  tube  bearing  much  of  the  load. 
In-plane  deformation  of  the  ring  cross-section  in  this  region  is  mostly  ovaling  (n  “ 2) 
with  some  n ” 3 mode  contribution,  for  the  short  time  considered  before  arrival  of  the 
overburden  load  from  atop. 
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7.2.2. 


Oblique  Front-On  Loading 


In  the  elastic  calculation  the  deformation  mode  in  the  yaw  plane  shows  the  longi- 
tudinal compression  accompanied  by  displacement  of  the  headworks  with  the  blast.  Bending 
of  the  headworks  in  the  yaw  plane  is  first  toward  the  blast  because  of  the  arrival  time 
effect  of  the  load  at  front  face  and  is  then  reversed  due  to  the  effect  of  loads  acting 
at  the  blast  side  of  the  headworks  which  is  exposed  for  this  calculation.  The  yaw  bending 
moment  induced  in  the  tubular  section  is  about  ono-half  of  the  pitch  piano  moment  for  the 
distribution  of  airblast  loads  assumed.  The  shear  stress  in  cross-sections  of  the  head- 
works  perpendicular  to  the  longitudinal  axis  does  not  seem  to  be  greatly  affected  by  the 
additional  out-of-plane  load. 

7.2.3.  Side-On  Loading 

The  deformation  mode  due  to  side-on  loading  under  the  two-dimensional  assump- 
tion wherein  only  the  tubular  section  of  the  shelter  is  analyzed  as  a ring  consists  mainly 
of  ovallng  with  some  breathing  mode  contribution. 

7.3.  AIRBLAST  CHARACTERISTICS  AND  DEFORMATION  MODES 

The  conclusions  of  airblast  characteristics  are  based  on  limited  airblast  data 
and  mostly  flastic  results. 

7.3.1.  Spatial  Characteristics 

Because  of  the  direct  application  of  load  to  the  exposed  surfaces  in  the  shelter 
concept  the  airhlast  overpressure  distribution  about  the  exposed  surface  has  significant 
effects  on  the  resulting  structural  response.  This  is  especially  true  because  of  the  com- 
plex geometry  involved  and  the  sensitivity  of  the  overpressure  characteristics  to  this 
geometry. 

In  the  front-on  loading  case  the  single  most  important  parameter  is  of  course 
the  reflection  coefficient  related  to  the  front  face.  Peak  stresses  in  the  headworks  and 
tube  are  directly  governed  by  the  peak  front  load.  The  load  distribution  about  the  front 
face  is  another  significant  factor  which  affects  the  amount  of  rotation  Imparted  to  the 
headworks.  Its  significance  is  diminished  for  an  inclined  front  face.  The  main  action 
of  the  overburden  load  which  is  to  push  the  headworks  downward  is  mostly  indirect  in  that 
the  overpressure  wave  has  to  transverse  the  backfill  medium  to  reach  the  structure  and 
does  not  exhibit  as  much  sensitivity  as  the  front  load.  The  magnitude  of  the  overburden 
load  is  also  substantially  less  than  the  front  load,  making  it  of  secondary  importance. 
Loading  on  the  sides  of  the  headworks  contributes  to  the  confinement  pressure  and  its 
effect  is  probably  significant  only  in  the  oblique  incidence  case  when  it  contributes  to 
deformation  of  the  structure  In  the  yaw  plane. 
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Another  front-face  load  parameter  which  merits  consideration  is  the  pressure  dis- 
tribution in  the  gap.  Its  magnitude  governs  shear  deformation  in  the  inner  upper  and  lower 
corners  of  the  frame  and  in  addition  provides  confinement  to  the  closure  concrete.  The  ef- 
fect of  its  distribution  in  the  gap  space  is  not  studied  here.  It  appears  that  the  effect 
of  minor  deviations  from  uniform  distribution  will  not  be  detected  in  the  response.  In  any 
case  the  impact  of  the  gap  on  the  frame  and  closure  is  localized  to  these  regions  and  is 
not  felt  in  the  headworks. 

In  contrast  to  the  direct  load  the  effect  of  indirect  loads  such  as  those  acting 
on  the  berm  is  moderated  by  passing  through  the  dissipative  backfill  material.  Any  varia- 
tion in  its  spatial  distribution  is  tempered  by  material  damping  by  the  time  the  ground 
shock  reaches  the  structure.  This  insensitivitv  is  indicated  in  the  two-dimensional  cal- 
culation where  a two-fold  increase  in  the  load  on  the  blast  side  of  the  berm  does  not  cause 
any  significant  increase  in  the  structural  response  as  measured  by  peak  bending  moment  and 
hoop  stress. 

7.3.2.  Temporal  Characteristics 

Among  the  many  characteristics  which  fall  under  this  category  the  most  important 
one  is  probably  the  impulse  associated  with  the  front  load  for  the  front-on  loading  condi- 
tion. When  the  closure  deforms  In  the  punching  mode  the  impulse-time  history  becomes  the 
parameter  which  governs  the  closure  response.  The  magnitude  of  the  total  impulse  and  the 
rate  at  which  this  impulse  is  applied  have  a direct  impact  on  the  maximum  and  permanent 
deformation  of  the  closure.  This  is  in  contrast  to  the  situation  when  the  closure  remains 
elastic  where  the  important  closure  response  considerations  are  vibration  and  resonance 
and  the  Important  airblast  characteristics  the  frequency  content  at  the  natural  frequency 
of  the  closure. 

In  general,  the  remark  made  earlier  about  spatial  characteristics  can  also  be 
made  about  the  temporal  characteristics,  l.e.,  when  the  loading  is  direct  the  transfer  of 
the  temporal  characteristics  of  the  load  to  structure  is  most  complete,  especially  the 
high  frequency  components.  However,  when  loading  is  not  direct  the  backfill  and  soil  mass 
surrounding  the  structure  will  act  to  filter  out  any  high  frequency  content  of  the  over- 
pressure and  only  the  lower  frequency  components  are  expected  to  be  important. 

Arrival  time  or  rather  the  sequencing  of  arrival  times  of  the  overpressure  at 
different  points  of  the  structure  is  also  another  important  characteristic.  In  the  front- 
on situation  phasing  of  the  front  load  and  overburden  load  determines  the  a.aount  of  rota- 
tion Imparted  to  the  headworks.  Arrival  time  differential  in  the  fr  t loads  in  the 
oblique  Incidence  case  determines  the  initial  bending  moment  into  the  blast.  These  effects 
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are  transmitted  down  the  shelter  to  the  transition  section  and  tube  section  where  they  be- 
come really  important. 

7.4.  EFFECT  OF  INELASTICITY 

As  far  as  response  of  the  shelter  is  concerned  the  foremost  manifestation  of  the 
effect  of  inelasticity  is  in  the  closure  deformations,  l.e.,  the  punching  phenomenon  dis- 
cussed in  Section  5.1.  This  response  is  sensitive  to  the  strength  of  the  concrete  assumed. 

Other  inelasticity  effects  are  observed  in  the  frame  notably  the  inner  corners, 
in  the  bearing  ring  and  concrete  region  behind  it  in  the  headworks,  and  in  the  fore  portion 
of  the  tubular  section  near  the  headworks.  These  effects  are  mostly  confined  locally  to 
the  region  mentioned  and  do  not  influence  response  in  other  parts  of  the  shelter.  The  only 
exception  is  the  inelastic  deformation  in  the  closure  bearing  region  whereby  the  headworks 
deformation  is  uncoupled  from  the  effect  of  front  load  acting  at  the  closure. 

Limiting  the  strength  of  the  concrete  in  the  tubular  section  also  limits  the 
amount  of  axial  load  the  tube  section  can  transmit  and  thus  facilitates  the  uniform  dis- 
tribution of  axial  compression  in  the  section.  The  impact  of  this  effect,  if  any,  has  not 
been  identified. 

Inelasticity  in  the  medium  has  the  effect  of  modifying  the  stress  pulse  before 
it  reaches  the  structure  and  provides  material  damping.  It  follows  then  that  its  effect 
is  proportional  to  the  expanse  of  medium  separating  the  structure  from  the  load.  Hence 
it  can  be  surmised  that  the  most  significant  effect  of  inelasticity  in  the  medium,  aside 
from  its  effect  on  the  overburden  load  described  earlier,  is  to  diminish  long-term  re- 
sponse and  affects  only  the  gross  response  of  the  shelter. 

7.4.1.  Elastic  Versus  Inelastic  Analysis 

Having  identified  the  effect  of  inelasticity  on  the  response  of  the  structure, 
the  merits  of  elastic  analyses  performed  for  the  same  structure  can  be  assessed. 

By  themselves  elastic  analyses  of  the  shelter  can  yield  useful  (and  exaggerated) 
Information  on  (a)  possible  deformation  modes  of  the  structure  under  different  loading 
conditions;  (b)  sensitivity  of  these  deformation  modes  to  variations  in  overpressure  dis- 
tribution and  wave  form;  (c)  possible  areas  of  distress  which  need  further  examination; 
and,  of  course  (d)  stress  states  in  parts  of  the  shelter  where  the  material  remains 
elastic. 

When  used  in  conjunction  with  engineering  judgment  the  quantitative  data  obtain- 
ed in  elastic  analyses  can  be  extended  to  yield  qualitative  data  on  inelastic  behavior. 

When  used  in  conjunction  with  inelastic  analyses  an  elastic  analysis  serves  as  an  excellent 
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baseline  to  be  referenced  and  compared  to  in  order  to  assess  and  Identify  certain  effects 
of  inelasticity  and  their  sensitivity. 

An  elastic  analysis,  of  course,  will  not  provide  data  pertaining  to  deformation 
and  mechanism  which  are  purely  inelastic  in  nature,  e.g.,  closure  response.  Nor  will  it 
provide  good  estimates  for  long-term  effects  due  to  the  lack  of  damping  and  it  invariably 
overestimates  peak  values  for  short-term  response. 

The  reason  an  elastic  analysis  is  more  useful  here  than  in  other  situations  is 
that  the  shelter  structure  is  semi-buried  and  subjected  to  direct  loading  which  dominates 
the  short-term  response.  It  would  be  less  effective,  for  instance,  if  the  response  of 
the  tubular  section  under  side-on  loading  is  sought. 

In  the  current  application  inelastic  analysis  also  has  its  limitations.  Fore- 
most is  the  uncertainty  of  material  properties  (modeling,  parameters,  etc.),  especially 
in  the  dynamic  realm.  Dynamic  plain  concrete  properties  are  not  readily  available,  not 
to  mention  reinforced  concrete  properties  which  reflect  the  complex  reinforcing  patterns 
in  the  shelter  configuration.  Soil  properties  have  always  been  a major  difficulty  in 
structure-medium  interaction  analyses,  although  in  the  case  of  the  shelter  the  interac- 
tion effect  is  not  critical. 

7.5.  GROUND  MOTION  EFFECT 

The  effect  of  airblast  induced  ground  motion  on  shelter  response  has  not  been 
studied  and  hence,  definitive  conclusions  on  this  aspect  of  the  dynamic  analysis  and 
the  role  of  the  soil  Island  approach  in  shelter-medium  interaction  calculations  cannot 
be  made. 

Because  direct  airblast  loading  dominates  the  shelter  response,  in  the  short 
term  at  least,  it  is  expected  that  the  ground  motion  effect,  if  any,  will  affect  only 
the  long-term  response  and  the  gross  behavior  of  the  shelter  structure  such  as  the 
longitudinal  bending  and  in-plane  deformation  of  the  tube  but  not,  say,  the  closure 
deformation. 

Based  on  ground  motion  results  obtained  in  Weidlinger  Associates'  study,  DNA 
001-77-C-0036,  there  is  no  appreciable  outrunning  ground  shock  in  the  (revised)  baseline 
site  at  the  600  psi  range.  In  the  outrunning  region  (100  psi  to  20  psi),  the  motions  ex- 
hibit a sustained  low  frequency  (1.5  Hz)  response  in  both  the  horizontal  and  vertical 
directions.  Ground  motion  effect,  therefore,  is  expected  to  influence  the  response  of  a 
shelter  located  at  this  range,  the  direct  airblast  effect  no  longer  as  dominant  due  to 
the  lower  overpressure  involved. 
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Depth  of  water  table  (100  to  300  feet)  is  found  not  to  have  significant  effect  on 
the  surface  ground  nK}tion.  The  effect  of  water  table  at  depths  less  than  100  feet  has  not 
been  studied. 

From  the  above  discussions  it  can  be  concluded  that  the  direct  ground  motion  ef- 
fect on  shelter  response  is  probably  secondary  to  the  airblast  effect  for  the  site  and  en- 
vironment considered.  However,  if  the  site  profile  is  changed,  for  example  by  raising  the 
water  table,  a reassessment  of  this  tentative  conclusion  will  be  necessary. 

From  a numerical  analysis  viewpoint  the  use  of  appropriate  ground  motion  at  the 
soil  island  boundaries  effectively  extends  the  size  of  the  Island  thereby  enabling  dynamic 
calculations  to  be  carried  out  for  a longer  period  than  would  be  possible  otherwise.  This 
is  absolutely  necessary  in  large-scale  three-dimensional  calculations  such  as  those  attempt- 
ed in  the  present  study  when  long-term  structural  response  of  the  shelter  is  analyzed. 

7.6.  SHELTER  ANALYSIS  METHODOLOGY 

From  the  foregoing  it  appears  that  a three-dimensional  dynamic  analysis  of  the 
shelter  using  numerical  methods  can  best  be  performed  in  two  parts:  short-term  and  long- 
term. In  the  short-term  analysis,  response  of  the  front  end  of  the  shelter  which  is  ex- 
posed to  direct  airblast  load  is  examined.  Peak  stresses  in  the  headworks,  frame  and 
closure  and  peak  deformation  of  the  closure  are  obtained.  Lonf-term  analysis  includes 
what  happens  afterwards  and  addresses  mostly  the  response  of  the  tube  and  the  gross  motion 
of  the  shelter.  (In  the  rear-on  incidence  case  the  short-term  response  would  include  the 
time  from  first  arrival  at  the  rear-end  to  arrival  at  the  front.) 

This  distinction  of  short  and  long-term  responses  follows  directly  from  the 
shelter  concept  wherein  the  structure  is  subjected  both  to  direct  airblast  and  indirect 
ground  motion  loads.  Detailed  subcomponents  of  the  shelter  are  located  up  front  (e.g., 
frame,  gap,  closure)  and  their  response  Is  dominated  by  the  direct  airblast  load.  Their 
important  response  time  Is  both  short,  about  the  order  of  magnitude  of  the  main  pressure 
pulse  (up  to  30  msec,  at  the  most  for  the  environment  considered  here),  and  local.  The 
fundamental  period  of  the  closure  for  Instance  is  about  8 msec.,  whereas  the  period  for 
longitudinal  bending  of  the  tube  is  100  msec.,  and  ovaling  of  the  ring  cross-section  about 
40  msec. 

The  distinction  can  also  be  made  from  a frequency  viewpoint.  High  frequency 
effects,  if  any,  are  more  likely  to  come  from  the  direct  airblast  load.  Airblast  induced 
ground  motion  by  its  nature  is  more  appropriately  called  low  frequency  input. 
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The  short  and  long-term  response  distinction  Is  also  one  of  (modeling)  analysis 
expedlancy.  Modeling  of  the  front  end  Is  by  necessity  a detailed  task  In  order  to  repro- 
duce the  response  characteristics  of  the  subcomponents,  say,  the  closure.  The  detailed 
model  usually  puts  a strain  on  the  capacity  of  the  computer.  The  core  storage  of  the  com- 
puter is  taxed  If  the  fine  mesh  used  for  the  front  end  is  continued  all  through  the  shelter/ 
berm  configuration.  Small  element  sizes  for  the  front  end  imply  small  Integration  time 
steps  for  an  explicit  code  (an  analysis  program  which  uses  an  explicit  solution  procedure) . 

Either  of  these  restraints  is  undesirable  and  is  to  be  avoided  if  possible. 

Results  obtained  In  the  present  study  indicate  that  the  two-part  analysis  metho- 
dology Is  both  practicable  and  correct.  The  key  observation  based  on  Phase  2 studies  Is 
that  minor  details  of  the  closure  do  not  significantly  affect  response  of  the  main  shelter 
and  vice  versa.  This  suggests  that  analysis  of  the  closure  and  main  shelter  can  be  un- 
coupled to  certain  extent,  leading  to  the  short-term  and  long-term  methodology  proposed 
here. 

An  analysis  procedure  for  performing  dynamic  three-dimensional  finite  element 
analyses  of  a shelter-like  structure  can  be  formulated  based  on  this  methodology.  The 
procedure  is  outlined  in  the  following. 

As  in  any  finite  element  analysis,  the  procedure  starts  with  modeling  considera- 
tions then  moves  on  to  guidelines  on  performing  the  dynamic  analysis  in  two  parts,  short- 
term and  long-term.  An  attempt  has  been  made  to  keep  the  comments  and  references  general 
enough  for  the  class  of  shelter-like  structure  considered;  specific  examples  are  made  to 

I i 

illustrate  a point  even  though  the  example  may  not  be  valid  for  a particular  shelter  con-  t 

figuration. 

I 1 

A . Modeling  Consideration  j 

Special  attention  in  modeling  should  be  given  to  the  following  shelter  subregions: 

• Closure — sufficient  details  should  be  included  to  reproduce  the  basic  deforma-  t 

tion  mode  which  is  punching  and  shear  in  the  Inelastic  range  and  fundamental  period  in  the 

elastic  range.  Any  steel  back  and  side  plates  should  be  so  modeled  as  to  reflect  its  bend- 
ing capability.  Concrete  material  properties  should  be  as  accurate  as  possible  since  the  ^ 

closure  response  (but  not  the  main  shelter  structure)  is  sensitive  to  the  material  properties  f 

used.  ! 

• Front  Face — the  front  face  model  should  be  as  detailed  as  the  air  overpressure  j 

distribution  on  it  dictates  in  order  to  bring  out  any  axial  load  imbalance  which  acts  to  J 

rotate  the  headworks.  1 
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• Frame/Gap — existence  of  the  gap  (and  gap  pressure)  must  be  Included  in  the  model. 
The  inner  lower  and  upper  corners  of  the  frame,  for  a square  closure,  as  well  as  any  other 
steel  lining  on  the  frame  must  be  modeled. 

• Bearing  Ring — a high  stress  gradient  exists  in  the  bearing  area  supporting  the 
closure.  A fine  mesh  is  necessary  here  to  reproduce  the  stress  concentration  that  is  the 
bearing  load.  Yielding  and  deformation  of  the  bearing  region  will  affect  the  overall  re- 
sponse of  the  closure. 

• Transition  Section — the  high  axial  stress  from  the  front  load  is  amplified  as 
the  load  bearing  area  decreases  and  a high  stress  gradient  results  in  this  region.  A moder- 
ately detailed  model  for  the  transition  section  is  probably  sufficient  although  any  abrupt 
change  in  cross-sectional  shape  (such  as  those  in  SI  and  SA)  should  be  modeled  with  care. 

• Fore  Tube  Section — the  front  portion  of  the  tubular  section  immediately  behind 
the  headworks  bears  the  axial  compression  and  bending  moment  transferred  from  the  headworks. 

A fine  subdivision  of  elements  should  be  used  in  this  critical  region  and  at  the  crown  in 
particular.  It  should  also  be  sufficiently  detailed  to  accommodate  in-plane  deformation 
modes  as  high  as  n = 3. 

• Main  Tube  Section — the  model  for  the  main  tubular  section  of  the  shelter  should 
reflect  its  dynamic  characteristics  such  as  longitudinal  flexure  and  in-plane  bending. 

• Soil  Medium — careful  modeling  of  the  soil  medium  immediately  in  front  of  the 
shelter  is  necessary  not  so  much  to  admit  correct  transmission  of  the  airblast  load  as  to 
prevent  erroneous  distribution  of  this  load  to  the  structure.  Modeling  of  soil  medium  else- 
where should  be  consistent  with  the  dynamic  characteristics  of  the  tube  section. 

B.  Airblast  Characteristics 

Special  care  should  be  exercised  in  determining  the  following  airblast  character- 
istics since  they  have  significant  influence  on  structural  response.  The  exact  effect,  of 
course,  varies  with  loading  conditions.  In  general  the  characteristics  are  more  important 
in  the  event  of  a direct  loading  than  otherwise. 

• Head-On  Condition — airblast  parameters  of  Importance  are  the  front-face  reflec- 
tion factor,  its  spatial  distribution  about  the  front-face  and  Impulse  and  Impulse  time 
constant.  The  reflection  coefficient  governs  the  peak  pressure  and  hence  the  peak  stresses 
in  the  headworks  and  closure  Its  distribution  about  the  front  face  is  also  Important  in 
determining  the  amount  of  bending  of  the  headworks.  Impulse  is  Important  to  the  closure 
when  the  latter  deforms  in  the  punching  mode  and  the  peak  inelastic  deflection  is  then 
determined  by  the  Impulse. 
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• Oblique  Front-On  and  Slde-On  Condition — reflection  coefficient  at  the  blastward 
and  rear  sides  of  the  headworks,  if  exposed,  is  an  important  parameter.  It  governs  the 
amount  of  bending  in  the  yaw  plane  and  its  importance  increases  as  it  approaches  the  front- 
face  reflection  coefficient  in  value. 

The  effect  of  pressure  distribution  on  the  gap,  especially  the  differential  be- 
tween the  blastward  side  and  the  rear  sides,  needs  to  be  assessed. 

• Overburden  Load — effect  of  overburden  load  on  the  shelter,  whether  it  is  head- 
on,  side-on  or  rear-on,  is  not  expected  to  be  as  significant  as  the  direct  load  parameters. 
Its  action  is  to  push  the  structure  down  and  with  the  blast  and  is  usually  tempered  some- 
what by  the  expanse  of  backfill  surrounding  the  structure.  Its  characteristics  should  be 
as  refined  as  the  backfill  properties  and  tube  response  characteristics  require. 

C.  Short-Term  Analysis 

A short-term  analysis  of  the  response  of  the  front-end  should  be  based  on  a fine 
model  of  the  closure,  front  face,  frame/gap,  bearing  ring  and  transition  section  (see  Sec- 
tion 7.1).  It  may  be  desirable  to  begin  with  an  elastic  analysis  which  can  serve  as  a 
baseline  for  the  subsequent  inelastic  calculations. 

That  portion  of  the  shelter  after  the  transition  section  may  be  replaced  by  a 
simpler  model  which  can  reproduce  approximately  (i)  the  impedance  to  axial  load;  (ii)  the 
energy  radiation  along  the  longitudinal  axis;  and  (ill)  the  resistance  to  bending  of  the 
headworks.  Springs  and  dashpots  for  Instance  can  be  used  to  replace  the  aft-half  of  the 
tubular  section  and  medium  in  the  longitudinal  direction  while  holding  it  fixed  in  the 
transverse  direction.  The  short-term  analysis  should  be  terminated  after  peak  stresses 
in  the  headwork  and  frame/closure  and  maximum  deflection  in  the  closure  have  been  obtained. 

An  inelastic  model  of  the  closure  should  be  used  and  it  should  be  sufficiently 
refined  to  produce  the  punching  deformation  mode  and  the  shear  crack  propagation  phenomenon 
for  sufficiently  high  loading. 

D.  Long-Term  Analysis 

Long-term  analysis  includes  also  the  structure  response  from  wave  front  arrival; 
but  it  is  not  designed  to  yield  response  of  detailed  substructures  (e.g.,  closure)  pro- 
vided by  the  short-term  analysis.  It  should  therefore  be  based  on  a course  model  of  the 
front  face,  closure,  frame/gap  and  bearing  ring,  but  using  the  appropriate  models  for  the 
transition  section,  tube  section  and  medium  as  described  in  Section  7.1.  Care,  however, 
must  be  used  in  choosing  the  approximate  front  face  model  so  that  Important  overpressure 
characteristics  are  not  lost.  The  approximate  closure  model  also  must  be  choosen  so  that 
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it  Is  still  capable  of  transferring  the  front  face  load  to  the  headworks.  For  instance, 
the  composite  closure  of  the  S4  configuration  may  be  replaced  by  an  equivalent  all  reinforc- 
ed concrete  model  for  analysis  purposes  so  that  the  steel  back  and  side  plates  no  longer 
Impose  constraints  on  the  time  step  requirement. 

To  carry  the  calculation  out  for  a long  period  (about  200  msec,  for  the  S4  case) 
the  soil  island  approach  must  be  used  to  limit  the  Interference  due  to  unwanted  signals  com- 
ing from  the  boundaries.  The  extent  of  error  Introduced  by  boundary  signals  depends  on  the 
loading  condition,  size  of  the  soil  Island  and  dissipative  properties  of  the  medium.  An 
elastic  analysis  is  probably  not  appropriate  for  a long-term  analysis  mainly  because  it 
lacks  material  damping  characteristics.  Ground  motion  input  at  the  soil  Island  boundaries 
is  definitely  necessary  in  this  case.  In  general  an  elastic  analysis  tends  to  overestimate 
the  response  and  Introduce  spurious  signals. 
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APPENDIX  A 


MODIFICATIONS  TO  TRANAL 

TRANAL  is  designed  to  treat  three-dimensional  soil-structure  systems  subjected  to 
airblast  and  induced  ground  shock  loadings.  The  structure  is  assumed  totally  burled  and  the 
alrblast  loading  surface  a flat  ground.  The  air  overpressure  is  fashioned  after  Brode  and 
is  retained  in  its  analytical  form  in  the  family  of  subroutines  BR73. 

In  preparation  for  the  planned  soil-structure  investigation  and  the  S4  pretest 
calculation  extensive  modifications  have  been  made  to  the  code  to  facilitate  the  rather  com- 
plex geometry  and  applied  load  distribution.  The  airblast  loading  surface  is  no  longer  flat 
nor  a single  surface  and  the  overpressure  time  histories  are  given  in  digitized  format.  In 
addition,  because  of  the  size  of  the  problem  Involved,  other  modifications  to  the  code  are 
also  necessary  in  order  to  fit  the  maximum  size  into  the  computer.  These  changes  will  be 
described  briefly  in  this  appendix  which,  however,  is  not  meant  to  be  a formal  documentation 
of  the  code.  Most  of  the  necessary  programming  is  done  by  John  Baylor  of  Weldllnger  Associ- 
ates, New  York. 

A.l.  APPLIED  PRESSURE  LOAD  MODIFICATIONS 

Pressure  load  can  be  applied  to  the  exposed  surfaces  of  the  structure  and  berm  in 
one  of  two  ways  depending  on  the  details  of  the  geometry  and  overpressure  involved  and  of 
course  the  users'  choice. 

A.  1.1.  Overall  Loading 

This  method  of  applied  load  is  Intended  for  large  surface  areas  (such  as  the  berm 
surfaces  or  any  large  expanse  of  exposed  surface) . Where  the  geometry  is  simple  (plane)  and 
the  pressure  load  does  not  change  drastically  from  one  location  to  another.  This  loading 
mode  constitutes  a generalization  of  the  code's  original  applied  pressure  (Brode)  module. 
Discretized  pressure  time  histories  are  given  for  points  on  a spatial  network  covering  the 
surface  to  be  loaded.  The  required  loading  at  any  point  on  the  surface  is  obtained  through 
spatial  interpolation  of  the  data  Involved. 

A. 1.2.  Detail  Loading 

This  mode  of  applied  load  is  tailored  to  complicated  surface  description  (such  as 
those  on  the  front  face,  headworks  and  gap  in  the  structure,  neighboring  soil  surfaces  and 
cutouts  in  the  driveway)  and  Irregular  pressure  distribution.  In  its  most  refined  form, 
each  discretized  structural  element  can  have  its  own  pressure  loading  description  (or  pres- 
sure cell)  consisting  of  time-pressure  pairs,  loading  direction  and  time  of  arrival. 
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A. 1.3.  Code  Modifications 


Changes  to  the  code  are  made  to  the  input  module,  element  table  set  up  and  applied 
load  processing.  The  modified  and  new  subroutines  are  described  briefly  below  approximately 
in  their  order  of  execution;  names  with  an  asterisk  correspond  to  new  subroutines. 

READ  PROCESSOR 

UNT  OPT  — expanded  to  include  input  data  unit  option  for  pressure  cell  specifica- 
tion 

PRT  OPT  — expanded  to  Include  option  to  print  pressure  cell  data 

CAP  75R  — material  property  specification  is  expanded  to  include  pressure  cell 
"property,"  with  zero  mass  and  scale  factor  and  arrival  time  substi- 
tuted for  elastic  bulk  and  shear  moduli 

CHK  LCL*  - pressure  time  history  number  is  checked  for  validity 

RED  LOD*  - pressure  data  are  read  in 

RED  TBL*  - pressure  time  history  pains  are  read  in 

CHK  LOD*  - pressure  data  is  checked  for  validity 

ELM  SET  — new  name  for  old  block  of  codes  in  SRTOUT 

ELM  MAT*  - pressure  cell  subroutine 

ELM  RED*  - pressure  cell  specification  data  are  read  in 

CHK  ELM  — cell  specification  is  checked  for  validity 

CHK  MAT*  - pressure  cell  specification  is  checked  for  validity 

CNG  NAM*  - utility  routine,  name  of  a data  block  is  changed 

EUi  SRT*  - utility  routine,  elements  are  rearranged  according  to  processing 
(zone)  order 

BLD  TBL  — element  table  set  up  routine,  modified  for  pressure  cell  applications 

ILZ  ZON*  - utility  routine,  sets  up  zone  count  parameters 

ILZ  ELM*  - utility  routine,  sets  up  element  count  parameters 

SET  ELM*  - element  table  entry  is  updated  for  pressure  cell  applications 

MAT  PTR  — material  pointer  is  inserted  in  element  table 

RUB  MAT*  - utility  routine,  a data  block  is  erased 
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COORDINATE  PROCESSOR 


B MATRIX  - slight  modification  to  accommodate  pressure  cell  specification 

NOD  MAS  — slight  modification  to  accommodate  pressure  cell  specification 

XEQ  PTR  — expanded  to  set  up  pressure  cell  related  parameters  for  execution 
phase 

TIME  PROCESSOR 

ZP  HEX*  — expanded  to  include  new  pressure  options 

CHK  LCE*  - element  to  be  processed  is  checked  to  see  if  it  is  a pressure  cell 

UPL  CNF*  - nodal  forces  for  a pressure  cell  are  updated 

LOD  CELP*  - pressure  in  a pressure  cell  is  computed 

LOD  CELF  - forces  in  a pressure  cell  is  computed 

UP  LCNV*  - nodal  velocities  for  a pressure  cell  are  updated 

STRS  BD  — modified  for  generalized  overpressure  load 

MAS  CHK*  - mass  of  boundary  node  is  checked  for  non-zero  value 

FC  CNTR*  - utility  routine,  determines  centroid  of  a surface  element 

FELIX  W*  - generalized  overpressure  load  is  determined 

WONG  FX*  - generalized  overpressure  load  is  determined 


CHK  IDX*  - utility  routine,  surface  index  for  pressure  cell  is  checked  for 
validity 

DMP  PAK*  - utility  routine,  a packed-word  data  block  is  unpacked  and  dumped 

NORMAL*  — utility  routine,  unit  normal  vector  to  a surface  is  determined 

PAK  WRD*  - the  six  surface  itidices  are  packed  and  meshed  into  the  element  table 
entry  for  each  element 

UNP  WRD*  - unpacked  the  word  packed  in  PAK  WRD 

A pressure  cell  is  designated  by  an  element  with  zero  mass  and  has  a corresponding 
material  description  with  material  number  equal  to  the  pressure-time  history  number  to  be 
applied.  As  such  it  is  equivalent  to  specifying  the  material  property  of  an  element  and  is 
Indeed  used  for  that  purpose.  Since  ordinarily  the  material  description  in  TRANAL  is  des- 
cribed in  terms  of  region  (a  region  can  have  more  than  one  element,  but  all  elements  within 
a region  must  have  the  same  material  property).  This  affords  the  user  another  option  in 


specifying  material  makeup  of  a region  and  has  been  found  to  be  quite  useful  especially  when 
only  a few  elements  of  an  existing  model  need  to  be  changed. 

Prescribing  the  pressure  load  within  an  alement  is  only  part  of  the  pressure  cell 
description;  additional  information  must  be  given  to  specify  which  of  the  six  faces  which 
make  up  an  element  are  acted  upon  by  this  load.  This  is  done  by  specifying  any  or  all  of 
the  surface  indices  1,  2,  3,  4,  5 or  6.  These  data  are  then  packed  into  one  word  and  stored 
as  an  entry  to  the  element  table.  Before  each  subsequent  use  then  this  word  must  be  unpacked 
and  the  required  data  retrieved.  This  process  enables  the  original  structure  of  the  code, 
l.e.,  element  table  with  one  word  per  element,  be  conserved  although  a nominal  penalty  in 
execution  time  is  incurred. 

A.  1.4.  Sample  Check  Problem 

The  elastic  axlsymmetric  problem  of  a steel-lined  cavity  in  rock  is  used  as  a for- 
mal check  case  for  the  applied  pressure  load  modifications. 

Uniform  load  in  the  form  of  a triangular  pressure  pulse  is  applied  at  the  surface 
of  the  steel  liner.  Three  TRANAL  runs  are  made  using  the  model  of  the  steel  and  medium  as 
shown  in  Figure  A-1.  The  same  pressure  time  history  is  used  in  all  three  runs,  but  a differ- 
ent mechanization  is  exercised  each  time  corresponding  to  the  original  (old)  loading  module, 
the  overall  loading  module  and  the  detail  loading  module  (pressure  cell  method)  described 
previously. 

Results  obtained  from  the  three  runs  are  Identical.  The  TRANAL  code  results  have 
also  been  compared  with  numerical  results  obtained  by  Lockheed  Palo  Alto  Research  Laboratory 
and  this  comparison  is  given  in  Figure  A-2  for  membrane  stress  in  the  steel  liner  and  in 
Figure  A-3  for  the  radial  displacement  also  in  the  liner. 

Extensive  functional  checkout  of  the  routines  has  also  been  performed,  but  is  not 
described  herein. 

A. 2.  OTHER  MODIFICATIONS 

A.  2.1.  Region  Coordinate  Edge  Processing 

The  region  coordinate  edge  processing  is  changed  from  the  whole  soil  island  basis 
to  a zone  basis;  this  change  is  made  so  that  larger  soil  island  (problem)  can  be  accommodated 
in  the  machine,  at  least  through  the  COORDINATE  PROCESSOR  stage  and  that  the  approach  is  more 
consistent  with  the  subsequent  TIME  PROCESSOR,  which  also  processes  the  elements  on  a zone 
by  zone  basis.  Four  additional  peripheral  units  need  to  be  declared;  one  for  coordinate 
and  the  other  three  for  x-edges,  y-edges  and  z-edges,  respectively,  if  these  are  not  linear 
edges.  The  effected  subroutines  are; 
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READ  PROCESSOR 
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REG  CRD  — modified  to  accommodate  new  coordinate  edge  processing  approach 

REG  EDG  — modified  "o  accommodate  new  coordinate  edge  processing  approach 

COORDINATE  PROCESSOR 

COR  DPP  — modified  to  accommodate  new  coordinate  edge  processing  approach 

RGN  CRD  — modified  to  accommodate  new  coordinate  edge  processing  approach 

BLD  RGN*  - region  coordinates  for  a zone  are  built  and  saved 
BLD  RDG*  - region  edge  coordinates  are  built  and  saved 
A.  2.2.  Stability  and  Volume  Check 

For  an  explicit  code  using  the  finite  differencing  scheme  to  step  forward  in 
time  the  integration  time  step  used  must  satisfy  the  Courant's  stability  criteria,  l.e., 
the  time  step  used  must  be  a fraction  of  the  smallest  transit  time  across  an  element.  This 
check  which  prior  to  this  time  is  the  responsibility  of  the  user  is  mechanized.  At  the  same 
time  and  a negative  volume  check  is  also  made  to  detect  errors  in  coordinate  data  which  may 
result  in  a negative  volume  for  the  element.  The  modified  and  new  subroutines  involved  are: 

COORDINATE  PROCESSOR 

COR  DPP  — modified  for  stability  and  volume  check 
STB  ILZ*  - stability  check  zone  parameters  are  set  up 
DST  RCP*  - utility  routine,  determines  distance  between  points 
STB  CHK*  - the  stability  characteristics  of  an  element  are  checked 
STB  DMP*  - utility  routine,  the  stability  table  is  dumped 
STB  SET*  - parameters  for  stability  check  are  set  up 
A. 2. 3.  Region  Material  Description 

Prior  to  this  change  the  region  material  data  are  read  in  according  to  zone;  this 
has  the  disadvantage  that  any  change  to  the  zona  makeup  will  necessitate  a corresponding 
change  to  the  region  material  input  data.  This  modification  enables  the  material  data  to 
be  read  on  the  soil  Island  basis,  hence  Independent  of  zone  configuration.  The  modified 
subroutine  is  ZND  CRP.  Naturally  it  is  no  longer  possible  to  designate  a whole  zone  as 
VOID  and  the  necessary  VOID  regions  must  be  included  in  the  input  data. 
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